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Abstract 

We consider a heat problem with discontinuous diffusion coefficients and dis¬ 
continuous transmission boundary conditions with a resistance coefficient. For all 
compact (e, 5) -domains P C M n with a d -set boundary (for instance, a self-similar 
fractal), we find the first term of the small-time asymptotic expansion of the heat 
content in the complement of P, and also the second-order term in the case of a 
regular boundary. The asymptotic expansion is different for the cases of finite and 
infinite resistance of the boundary. The derived formulas relate the heat content to 
the volume of the interior Minkowski sausage and present a mathematical justifica¬ 
tion to the de Gennes’ approach. The accuracy of the analytical results is illustrated 
by solving the heat problem on prefractal domains by a finite elements method. 

Keywords: heat content; discontinuous transmission condition; Minkowski sausage. 


1 Introduction 

We consider a compact domain P C M n with boundary <9P that splits R" into “hot” 
and “cold” media, P + = P and P_ = R n \ P, characterized by (distinct) heat diffusion 
coefficients D + and D- (Fig. Q]). On the boundary f)P is also defined a function 0 < 
A (a;) < oo which describes the resistivity to heat exchange through the boundary. 

We are interested in propagation of the heat content associated with the following 

*Laboratory Jacques Louis Lions, University of Paris 6, Pierre et Marie Curie, 4 place Jussieu, Paris, 
France, claude.bardos@gmail.com 

lLaboratoire de Physique de la Matiere Condensee, CNRS - Ecole Poly technique, Palaiseau, France, 
denis. grebenkov@poly technique. edu 

■^Laboratory Applied Mathematics and Systems, CentraleSupelec Paris, Grande Voie des Vignes, 
Chatenay-Malabry, France, anna.rozanova-pierrat@centralesupelec.fr 


1 



problem: 


d t u± — D±Au± = 0 x G fl±, t > 0, 

u + \t=o = 1 , U-\ t =o = 0 , 


D_ 


D. 


du _ 

dn 

du + 

dn 


an — A(x)(m_ — u + )|an, 


Ian 


= Zb 


du _ 
<9n 


Ian, 


( 1 ) 

( 2 ) 

(3) 

(4) 


where d/dn is the normal derivative directed outside the domain hi. 

A rigorous analysis of the problem (|T |l - (ff|) for irregular boundaries requires its vari¬ 
ational formulation in appropriate functional spaces (see Section ED- The variational 
problem is shown to have a unique weak solution with the desired trace properties on 
the boundary dVt (see Section E])- The variational problem is equivalent to the prob¬ 
lem 0-0 for a piecewise Lipschitz Oil according to the classical trace theorem. In 
turn, extensions of the trace theorem have to be used for fractal boundaries or, more 
precisely, d-sets (see Subsection E2D- 

Once a unique solution u± of the problem (E])-(I3]) is established, we study the asymp¬ 
totic expansion of the heat content as t —> 0 


N(t) = 



U-(x, t)dx 


Vol(O) — / u + (x, t)dx. 


(5) 



Figure 1: Illustration of the heat content problem for a planar domain Q with prefractal 
boundary dfl presented by the third generation of the Minkowski fractal (of fractal 
dimension 3/2). This boundary splits the plane into two complementary regions. At 
time t = 0, the inner region h2 + = Q is “hot” (functions on h2 + are denoted with 
subscript + ), while the outer region f2_ = R n \ is “cold” (functions on are denoted 
with subscript —). 

Eqs. m-m describe heat exchange between two media prepared initially at differ¬ 
ent temperatures and separated by a partially isolating boundary[Tl [2], In fact, u(x,t) 
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can describe how the distribution of (normalized) temperature evolves with time. The 
transmission boundary conditions (151) . (151) impose the continuity of the temperature flux 
across the boundary, and relate this flux to the temperature drop at the boundary due 
to thermal isolation. The growth rate of the heat content with time characterizes the 
efficiency of thermal isolation. Understanding this problem is relevant to improve heat 
exchangers, e.g., cooling of metallic radiators or thermal isolation of pipes and buildings. 
Depending on application, cooling rate has to be either enhanced (e.g., in the case of 
microprocessors or nuclear reactors), or slowed down (e.g., in the case of pipes and build¬ 
ings). For these purposes, one can either modify the thermal isolation (i.e., the resistivity 
A), or the shape of the exchange boundary. It is therefore crucial to understand how the 
shape of the boundary influences heat exchange. In particular, would an irregular (e.g., 
fractal) boundary with a very large exchange area significantly speed up cooling? 

Similar equations can describe molecular diffusion between two media across semi- 
permeable membranes M- In that case, u(x,t ) represents the (normalized) concen¬ 
tration of molecules, while Eqs. o-o can model the leakage of molecules from a cell 
(h2 + ) to the extracellular space () or, more generally, the diffusive exchange between 
two compartments (e.g., oxygen or carbon dioxide exchange between air and blood across 
the alveolar membrane in the lungs). The resistance A is related to the cellular mem¬ 
brane permeability. As for heat exchange, one may need to enhance or to slow down the 
molecular leakage, and the shape of the boundary may play an important role. 

The discontinuity of the initial condition, of the diffusion coefficient, and of the solution 
u(x, t ) across the boundary between two domains constitutes one of the mathematical dif¬ 
ficulties to be treated. From a physical point of view, such discontinuities might appear 
unrealistic. For instance, the diffusive flux at the boundary at time t = 0 is infinite. 
For any physical setting of heat or molecular diffusion, there would be an intermediate 
layer between two media in which the material properties would change rapidly but con¬ 
tinuously. When the thickness of this intermediate layer is much smaller than the size of 
the domain, the physical problem with continuously varying parameters can be approx¬ 
imated by the heat problem (H])-((5j). Such an approximation is applicable starting from 
a small cut-off time while understanding the heat exchange at smaller time scales would 
need either restituting an intermediate layer, or introducing nonlinear terms into the heat 
equation. Throughout this paper, we focus on the mathematical problem ([T|)—(J4]). 

The physical properties of the two media and are supposed to be different: 
D + ^ D_ . This implies the discontinuity of the metric on dfl . The case of continuous 
metric (g_\gn = g+|an) on smooth compact n -dimensional Riemannian manifolds with 
a smooth boundary dfl was considered in Ref. |5j. The case of continuous transmission 
boundary conditions for the expansion of the heat kernel on the diagonal was treated 
in Ref. |6] (see also Ref. for a survey of results on asymptotic expansion of the heat 
kernel for different boundary conditions). The heat content asymptotic expansion with 
Dirichlet boundary condition was found 

• up to the third-order term for a compact connected domain D C M" with a regular 
boundary e C 3 (Refs. |8| [9]) ; 

• up to an exponentially small error for a compact connected domain fl C R 2 with a 
polygonal dQ (Ref. pTO] ) and for Ocl 2 with <9h2 given by the triadic Von Koch 
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snowflake (Ref. [HI]); 


• up to the second-order term for the general case of self-similar fractal compact 
connected domains in W 2 (Ref. [T2]). 

In general, the boundary between two media can have some resistance to heat ex¬ 
change, described by the function A(x) > 0 (x G <9f2) that may account for partial 
thermal isolation. We outline three cases of boundary conditions according to A ; 

1. If 0 < X(x) < oo for all x G <912, u is discontinuous on dfl and we have: 




an 


A(z)u+|an, 



an 


D. 


du- 

dn 


an- 


2. If A = +oo for all x G <912, u is continuous on <9f2 due to the transmission 
condition and in this case 


u + \an = u- |an, 



an 


D_ 


du _ 
dn 


an- 


3. If A = 0 for all x G <9Q, we have the Neumann boundary condition 

dii- du + 

-t :—\an = t:— an = 0 
dn dn 

that models the complete thermal isolation of <9f2 and implies the trivial solution 
given by U-(x,t) = 0 and u + (x,t ) = 1 for all time t > 0. 

The main goal of the article is to develop the preliminary study[T3] and especially 
to formalize the seminal approach by de Gennes[H]. In the case A = +oo, de Gennes 
argued that as t —> +0 , N(t) is proportional to the volume fi(dQ, y/D+t) of the interior 
Minkowski sausage of dVt of the width equal to the diffusion length y/D+t : 

fi(d£l,£) = Vol({a; G 01 dist(x,9I2) < £}) 

(see also Ref. mi). In particular, 

• for a regular boundary dfl, N(t ) is proportional to Vol(dfl)y/D + t ; 

• for a fractal boundary <912 of the Hausdorff dimension d, N(t) is proportional to 
( D+t)—. 

The de Gennes scaling argument was further investigated in Ref. [TB] , both experimentally 
and numerically. It was shown that irregularly shaped passive coolers rapidly dissipate 
at short times, but their efficiency decreases with time. The de Gennes scaling argument 
was shown to be only a large scale approximation, which is not sufficient to describe 
adequately the temperature distribution close to the irregular frontier. 

In the present paper, we provide a mathematical foundation and further understanding 
for the de Gennes approach. We obtain three results valid for all compact (e, 5) -domains 
12 in with connected boimdary <912, presented by a closed d-set (see Section [2~2l for 
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the definitions of (e, 5)-domains and d-sets): the well-posedness of the problem (1T]) - P]h 
the continuity of the solution on A (see Section [2]), and the asymptotic expansion of the 
heat content (0). In particular, these results hold for domains with a self-similar fractal 
boundary. 

We show in Theorem [5] that the heat content N ( t ) is approximated by the volume of 
the interior Minkowski sausage of of the radius y/AD + t: 


N(t ) = T\ 


C\(d£l)n (d£l, yj AD + t^j + O (dQ, yj AD+tj^ , 


( 6 ) 


where T\ is equal to 1 if A = oo and yft if A > 0 is finite. Here Cx(dQ) is a constant 
depending only on the shape of dfl and finiteness of A (see Theorem 0 for the exact 
formulas). Formula ([6]) is the first approximation of Eqs (j80|h (l82li given in Theorem 0 
which allows to find N(t) up to terms of the order T\0 (\/i // (<9fl, y/AD + t)) . 

Moreover, the asymptotic relation (J6]) remains valid even for mixed boundary condi¬ 
tions for three disjoint boundary parts, i.e. when A = oo on one part of the boundary, 
A = 0 on another part, and 0 < A < oo on the remaining boundary (see Theorem [31)- 
However, changes of the type of the boundary condition should be continuous (see The¬ 
orem ED such that u remains a continuous function of A. In this more general case, the 
coefficient C\ (<9fi) in Eq. (ED is given either by Eq. (l83li for 0 < A < oo, or by Eq. (1841) 
for A = oo, or is equal to 0 for A = 0 (the boundary with A = 0 does not contribute 
to the short-time asymptotics of the heat content). Finding the asymptotics for mixed 
boundary conditions with a discontinuous jump from a finite A to A = oo is still an open 
problem. 

As expected, the resistivity of the boundary to heat transfer makes heat diffusion 
slower due to the presence of the coefficient T\ = y/t. 

For a fractal boundary we replace /i (<9h2, y/AD + t) by the volume of the interior 
Minkowski sausage which scales as {AD + t)^ n ~ d ^ 2 , where d is the fractal dimension[T2], 
In the fractal case the integral over dfl should be understood by using the Hausdorff 
measure (see Ref. \m uni nrj). 

The comparison between the asymptotic formula (JHD and a numerical solution of the 
problem (0Q)-(]4D for the unit square and a prefractal domain is shown in Fig. [2] for a finite 
A and in Fig. [3] for A = +oo. The numerical solution was obtained in FreeFem-H§* 
by a finite elements method with the implicit 6 -schema, also known as Crank-Nicolson 
schema, for the time discretization with 9 = \ and At = 10~ 6 . The domain was 
centered in a ball B of diameter (at least) twice bigger than the diameter of H. The 
Neumann boundary condition was imposed on the boundary of the ball. According to the 
principle “not feeling the boundary” m (see also Section [3D, the heat content propagation 
in M 2 with a prescribed boundary dfl can be very accurately approximated at small times 
by the heat content propagation computed in B. The accuracy of this approximation 
can also be checked by changing the diameter of the ball. In the case of the square 
domain H, the ball was replaced by a square with four times bigger edge. Each pre¬ 
fractal edge was discretized with 27 space points while 57 points were used in the external 
boundary of the ball. The mesh size was varied to check the accuracy of the presented 
numerical solutions. For the case of the discontinuous solution on the boundary (when 
0 < A < oo ) we apply the domain decomposition method and match the boundary 


5 











Figure 2: Comparison between the asymptotic formula (J6]) (solid line) and a FreeFem f I 
numerical solution of the problem (JT])— (J4J) (circles) for two domains: (a) the unit square 
(Vol(<9f2) = 4) and (b) the third generation of the Minkowski fractal (Vol(<9f2) = 2 3 • 
4), with D + = 1/100, £)„ = 1, and A = 17. Since the Hausdorff dimension of the 
boundaries of these domains is 1 (even for the prefractal case), Eq. (J6]) for a constant 
A is reduced, according to Theorem 0 to N(t ) = 2\/tC 0 \fi(dQ, y/AD + t) + O(ti) with 
/i(<9Q, y/AD + t) ~ y/ADfl Vol(<90) and C 0 given by Eq. (1971) . For plot (b), dashed 
line shows the fractal asymptotic (that would be exact for the infinite generation of the 
fractal) with de Gennes approximation of /i (d£l, \JAD + t ) in Eq. (JHj) by (AD + t)i . This 
approximation is valid for intermediate times. 



Figure 3: Comparison between the asymptotic formula ([6]) (solid line) and a FreeFem—— 
numerical solution of the problem (lT]) - (14j) (circles) for two domains: (a) the unit square 
(Vol(<9f2) =4), and (b) the third generation of the Minkowski fractal (Vol(<9f2) = 2 3 • 4), 
with D + = 0.4, D_ = 1, and A = oo. Since the prefractal boundary has the 
Hausdorff dimension 1, Eq. (J5]) is reduced to Eq. (l95|) . i.e., N(t) oc \Jt. In turn, dashed 
line shows the fractal asymptotic (that would be exact for the infinite generation of the 
fractal) with de Gennes approximation of /i (<9f2, y/AD + t) in Eq. (EJ) by 2.5(4D + t)^ . This 
approximation is valid for intermediate times. 


values of the respective solutions on by a Picard fixed point method. We consider 
therefore the numerical solution of heat propagation for small times as a reference, to 
which asymptotic formulas are compared with. In particular, deviations between the 
numerical solution and the asymptotic formulas observed at longer times illustrate the 
range of validity of the short-time expansion. 
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For the regular case dfl G C 3 , we obtain the heat content approximation up to the 
third-order term. The formulas are given in Theorem [3[ For the case A < oo, the 

3 

coefficient in front of the second-order term (£2 ) in the asymptotic expansion depends 
on the mean curvature. In turn, for A = 00 , the second-order term (here, £) in the 
asymptotic expansion vanishes: 

m = 2l ~ e f7W^k Vo ' mVi+oith (7) 

The rest of the paper is organized as follows. I 11 Section [21 we describe the class of 
irregular boundaries and prove the well-posedness of the model relying on the variational 
formulation of the problem. The boundary conditions are treated in the weak sense by 
generalizing the trace operator and the Green formula to fractals using fractal Besov 
spaces, B^’ 2 (dtt) and B 2 j^{dVl) (/? = 1 — >0 for a d-dimensional dQ ) defined 

in 0 In Section [2] we also establish the continuity of u as a function of A . In Section [3] 
we prove that the problem to fold N(t) can be replaced by a heat problem localized 
in 0{\ft) -interior Minkowski sausage of the boundary by a variant of the principle “not 
feeling the boundary” m in the general case in R” . This allows, due to the continuity of 
it on A, to establish Theorem [3] for a mixed boundary condition including zero, finite, or 
infinite values of A . Considering a regular <9f£ (at least in C 3 ) and using the localization 
properties from Section [3J, we rewrite in Section [4] the formula for N (£) in the terms of 
the local coordinates. Section [5] gives the approximation of the heat problem solution 
through the solution of one-dimensional constant coefficient problem. The heat content 
is calculated in terms of the volume of the interior Minkowski sausage of the boundary in 
Section [31 Firstly, to illustrate the technique of the proof on a simple case, we give the 
proof for the case of continuous diffusion coefficients D + = ZA_ , just with discontinuity of 
the initial condition. I 11 this case, all formulas given in Section [ 6 ] are valid for all types of 
the boundary introduced in Subsection 12.21 The calculation relies on the Green function 
of the problem with constant coefficients for hi being a half-space (see [B]). We also obtain 
the Green function used in Section [7] for the proof of the asymptotic heat expansion up 
to the third-order term for a regular dVt G C 3 . 


2 Well-posedness of the model 

Let fl be an open connected bounded subset of R n such that dVt is closed with Vol(fi) < 
00 . We denote by = hi and = R n \ hi (Fig. [[]). 

We are looking for the solution of the problem (JT]) -(13]h where D + 7 ^ ZW , D + > 0 
and _D_ > 0, A(x) > 0 for all x G dfl . The boundary dfl is divided into two disjoint 
parts: Too = {x G <9fi| A(x) = + 00 } and dVL \ = {x G <9h2| 0 < A(x) < + 00 }. Each 

of the parts can be the empty set. We thus assume that A G L°°{dVt \ Too). 

2.1 Regular boundary: at least piecewise Lipschitz 

Firstly, we consider the case when dVt is regular (at least piecewise Lipschitz) and Too 
is the empty set. 


7 






To prove the existence, the uniqueness, and the stability of a solution of the prob¬ 
lem (HD—®, we P r °ceed with its variational formulation. 

We introduce the space H = L 2 (M. n ) and the space 

V = {feH\U = f\n + e H\n + ), and /_ = /| n _ (E 

of functions / = /+ ln+ + /-1 q_ defined on U such that their restrictions f + = 
f\n + and /_ = f\a_ belong to H 1 . We equip V with the norm: 

2 d x + D_ / |Vn_| 2 da;+ / |n| 2 da;. 

Jci- Jn + uo_ 

We notice that Id is a Hilbert space, V C L 2 (Q ), and V is dense in L 2 (Fl). In addition, 
V C L 2 (M n ) C V' , where V' is the dual space to V. Finally, since dfl is regular, the 
inclusion V C L 2 (R n ) is compact. 

Applying the usual trace theorem under the assumptions that is bounded and dfl 
is at least piecewise Lipschitz, the bilinear form 



a(u, v ) = D + 


S7u + S7v+ + D_ 


Vm_Vu_ 


'n_ 


is continuous, 



\(x)(u + 


uJ)(y + — v_)do (8) 


| a(u, n)| < C'(||A||Loo( a n), D + , D_, n + )||u||y||u||y (for a constant C > 0), 
and coercive on V x V , i.e., 


a(u,u) = D + / |Vu + | 2 dx + D_ / |Vn_; : da;+ / \(x)\u + — u_\ 2 da 


<dn + 


— 11 ^ 11 v 11 u 


> 0 . 


Thus we conclude(T 8 l| that the bilinear form a(u,v) defines an operator A : V —> V' 
by a(u,v) = (Au,v). Moreover, —A\ L 2 ^n^ with D(A) = {u E V\ Au G L 2 (R n )} 
generates an analytical semigroup. 


Remark 1 When Too is not empty, the variational form (G|) is well adaptable to the case 
where u is continuous across the part T^ C dil of the interface. By convention we put 
on this part X(x) = oo which implies u + = on T^ (see also Theorem^). 

For Too 7 ^ 0, we introduce V as the space of functions u G L 2 (M n ) such that 

u + — u \n+ G // 1 (fl + ), = u~\n_ G iF 1 (Q_), u + \ Foo = U-\ Foo , 

and, therefore, we consider the bilinear continuous and coercive form on V x V 


a{u,v) = D + / X7u + Vv + + D_ / Vu_Vu_ 


idn\ r 0 


A {x){u + —uJ){v + — vJ)da. (9) 


In particular, for Too = dfl, we get V = H 1 (X1 + U fi_) and 

a(u,v) = D + / Vu + Vv + -\-D_ / Vw_Vw_. 
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2.2 Extension to d -sets (fractal case for d > n — 1) 

Let us define a class of fractal domains to be considered. We will see that the existence and 
uniqueness results of a weak solution of the problem (TT]) -(pi]) hold for a class of bounded 
(e. h) -domains[20i 121. 22 ] f 2 + such that dLl is a d-set| 2 Tj: 

Definition 1 (d -set [21l 122 '. f 2Sf) Let r be a closed subset of M" and 0 < d < n. 
A positive Borel measure m<i with support T is called a d -measure of T if, for some 
positive constants C\, C 2 > 0 , 

c\r d < m d (T fl U r (x )) < c 2 r d , for V x G T, 0 < r < 1, 

where U r (x) C M n denotes the Euclidean ball centered at x and of radius r . 

The set T is a d -set if there exists a d -measure on T . 

As it is known from Ref. m, p.30], any two d -measures on I are equivalent. 

Definition 2 ( (e, 5) -domain l20i HO . \22hj ) An open connected subset Q of M” is an 
(e,S) -domain, e > 0 , 0 < 5 < 00 , if whenever x,y £ and \x — y\ < S, there is a 
rectifiable arc 7 C with length £( 7 ) joining x to y and satisfying 

1. £( 7 ) < and 

2. d(z, dLl) > e\x — for z £ 7 . 

I 11 particular, a Lipschitz domain is an (e, h) -domain and also a n-set|2 2 | | (i.e., a d-set 
with d = n). Self-similar fractals (e.g., von Koch’s snowflake domain) are examples of 
(e, 6) -domains with the d-set boundary |19l 22 ] . d > n — 1 . 

In order to describe irregular boundaries of fractal dimension d > n — 1, we define 
sets preserving Markov’s inequality (Ref. [2TJ Ch. II): 

Definition 3 A closed subset V in M” preserves Markov’s inequality if for every fixed 
positive integer k, there exists a constant c = c{V, n,k) > 0, such that 

c 

max VP < - max P 

vnu r (x) r vnu r (x) 

for all polynomials P G P& and all closed balls U r (x ) , x £ V and 0 < r < 1. 

Examples of sets that preserves Markov’s inequality are d -sets in W 1 , where d > n — 1, 
and self-similar sets that are not a subset of any (n — 1 ) -dimensional subspace of R n (see 
Refs. (221125]). 

To extend the variational formulation introduced in Subsection 12. II to fractal bound¬ 
aries of the type of d-sets, we use the existence of the d-dimensional Hausdorff measure 
m d on dfl (the d-measure from Definition [T]) and the theorem which generalizes the 
usual trace theorem and the Green formula. 

For example, for d = n — 1 and a Lipschitz dkl , we know(T 8 . m that the trace 
operator is linear continuous and surjective from H l {Lt) onto H^^dkl) , and the formula 

f du f 

J a vAudx = ~ j a VvVudx < 
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holds whatever u £ i7 1 (0) such that A u £ L 2 (fiV) and v £ i7 1 (0). 

To generalize the trace operator and the Green formula to fractal boundaries, one 
introduces the Besov space B 2 i ' 2 (dLl) with P = 1 - ^ > 0 (see Note that for 
d = n — 1 , one has (3 = \ and 

R?’ 2 (<9G) = if* (0ft), 

2 

i.e., one recovers the above relations. In general, 


1 . 


For an arbitrary open set ft of M n , the trace operator Tr is defined[2T1 [25 , [27] for 
u e L} oc {ti) by 


Tra(x) = lim ■ 


u{y)dy, 


( 10 ) 


r^om(nnu r (x)) J nnUr{x) 

where m denotes the Lebesgue measure. The trace operator Tr is considered for 
all x £ ft for which the limit exists. 


2. If ft is a bounded (e, 5)-domain in such that its boundary <9ft is a closed d-set 
preserving Markov’s inequality, then [211 [22 j 

a) the trace operator Tr : id 1 (ft) —* B 2 f 2 (dLl) is linear continuous and surjective; 


(b) the Green formula holds (see also Refs. [271 12S j for the von Koch 


case m 


,du 


J^vAudx = (^,Tru) ((B 2 , 2 (an))/B 2 , 2 (an)) - J^WvWudx, 

where the dual Besov space (B 2 f 2 (dLl))' = R 2 ’^(<9ft) is introduced in Ref. 
(see [Aj. 


( 11 ) 


Let us also notice that the Green’s formula (imi still holds whatever u £ id 1 (ft) such 
that Aw £ d7(ft) and v £ dd x (ft). 


2.3 Well-posedness 

The above preliminaries allow us to prove the following Proposition: 

Proposition 1 1 . Let ft = ft + be a bounded domain in M n with a closed piecewise 

Lipschitz boundary <9ft and 0 < A(x) < +oo be a given function defined on <9ft . 
By Too is denoted the part of <9ft such that 

Vx £ Too X(x) = +oo, 

in the such way that A £ dv°°(<9ft \ T^) . Then the bilinear form 

a(u,v) = D + / Vw + Vu + + D_ / Vw_Vu_ + / A(x)(w + — uJ){y + — v_)do 

JJ dfi \roo 

is continuous and coercive on V x V with 
V = {u £ L 2 {R n )| u + = u| n+ £ H\n + ), m_ = u\ Q _ £ 

u + = U- on Too}. ( 12 ) 
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2. Let be a bounded (e, 5) -domain in M n with a closed d-set boundary dLl 

and A G C(<9r2) be a positive continuous function defined on dLl. By is denoted 
the d -measure on dfl (see DefinitionUf). Then the bilinear form 


a(u,v) = D + / Vu+Vv-t 

J Q_|_ 


D_ / 


+ / A(x)Tr(u + — -u_)Tr(w + — v_)dmd 


i on 


is continuous and coercive on V x V (V is defined in Eq. H12\) ). 

3. Let LI = f2 + be a bounded (e, <5) -domain in M n with a closed d-set boundary dLl 
and A(x) = +oo for all x G dLl. Then the bilinear form 


a(u,v) = D + / Vit+Viq 

J Q_)_ 


D. 


Vw_Vu_ 




is continuous and coercive on V x V with V = H 1 {M n ) . 

Consequently, we obtain the following theorem: 

Theorem 1 (Well-posedness) In all cases from Proposition\]]for all u 0 G H = L 2 (W n ) 
there exists a unique solution u G L 2 (M”)) D L 2 (R)", V) of the variational problem 


Vn G V — (w, v )h + a(u, v ) = 0, u(x, 0) = uq G L 2 (R n ), 


(13) 


where by (■, •)# is denoted the inner product in H. In addition, this solution verifies the 
energy equality: 


1 

2 


9 f 1 

\u(t)\ 2 dx + / a(u,u)ds = - 

Jo 2 


\u 0 (x)\ 2 dx. 


(14) 


Remark 2 On one hand, any “smooth enough” solution of the problem w-w gives 
the solution of Theorem [2 On the other hand, any solution from Theorem . [7] satisfies 
the relations and, in a weak sense (in the sense of the duality presented above), 

satisfies the relations 

Finally, we prove 

Theorem 2 (Continuity of u\ on A and the case A = oo ) Let (A k)km be a positive 
sequence converging to A* in L°°(dLl ) . Then the corresponding sequence of the solutions 
(u\ k )keN of the system (|7])-Q) converges strongly to u* x in L 2 (R n )) fl l 2 (K t + , h), 

i.e., u\ is continuous as a function of A . 

If t oo in L°°(dLl), then u\ k —> u 00 in L 2 (R n ))flL 2 (R^, V) with (noo) + = 

(«oo)_ on dLl. In this case, E C(Rf, L 2 (W n )) H L 2 (Rf, solves 

\/v G i/ 1 (M n ) j d t u 00 vdx+ j D(x)'Vu 00 Vvdx = 0, u^x, 0) = uq G L 2 (R n ), (15) 

J R" Js. n 

with D(x) = tn+D+ + ln_D_ . 
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Proof. Firstly we suppose that A* is a finite bounded function on dfl (||A*||i / oo(gm < oo). 
Since u( 0 ) = u 0 does not depend on A, the equality (Thill implies that the sequence ( U \ k ) 
is bounded in C(Rf, L 2 (R n )) D L 2 (M t + , V). 

Therefore, due to the unicity of the solution for A* and the unicity of the weak limit, 
the convergence A*, —> X* in L°°(dLl) implies u\ k —*■ u \*. Since Uk\on G B 2 ’ 2 (dLl) with 
(3=1 — > 0 , with the help of (1T31) and the coercive behavior of a(u, u) , 

a(u,u ) > a||ti||y, for a > 0 , 

we conclude that u\ k —> u\* in C(Rf, L 2 (R n )) D L 2 (Rf, V). 

In the case || Afc|| ^oo^q) —>■ +oo, we find from (fT4l) that 

\/k eN / Afc(x)Tr((rtfc) + - ( u k )_) 2 dm d < oo, 

Jon 

{[ Tr (M+ - {u k )_) 2 dm^\ < 1 = [ \u 0 (x)\ 2 dx. 

\Jan J 

Therefore, we obtain in this case that (Woo)+ = { u oc)~ on 9Q, where by u Q0 we denote 
the limit of Uk as A —> +oo . In addition, u^t, •) G i/ 1 (M n ) and it is the solution of (1T5T) . 
□ 


3 Heat content localization to a small neighborhood of 
the boundary 


s the initial condition is zero in R n \ fl, we have 

N(t) = / (1 — u(x, t))dx = Vol(O) — f u(x,t)dx, (16) 

Jn Jn 

or equivalently, in terms of the Green function of the problem (JT ]) - (Ii]h 

N(t) = Vol(O) — / f G(x,y,t)dydx. 

Jn Jn 

Let us show that it is sufficient to integrate only on a small neighborhood of the boundary 
dLl to obtain the desired heat content with an exponentially small error: 


Lemma 1 Let F C Q be a non-empty open bounded set in M n ; such that dist(F, dLl) = 
e > 0 . Then for t —y +0 and u = u + \q + the solution of (Q])-(^]) ? associated 

with the Green function G(x, y, t ), 

1. it holds 



u + (x,t))dx = j 
= O 


1 -J G(x,y,t)dy J dx 

n—2 

p -e 2 /(4 D+t) 


VdD^t 


(17) 
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2. for e > 2 \fD+t such that e = 0(y/t), there exists 5 > 0 (a constant independent 
on time) such that the heat content N(t) can be expressed as 


N(t) — / U-(x,t)dx 
Jm. n \n 


1 — / G(x, y, t)dy ) dx + 0 [e <* ) , 

v «/ Qg 

where is the e -neighborhood of d£l. 


(18) 


Proof. As it was shown, the problem CQ)-(j4|) has a unique solution u = m+1q + 
. Let G(x,y,t) be the Green function so that 

u(x,t) = / G(x,y,t)dy. 

J Q 

Thus, using the properties of G such as G > 0 for all (x, y, t ) e M" x M n x M + and 
f Rn G(x, y, t)dy = 1, we easily see that 


0< [ G(x,y,t)dy = u{x,t) < [ G(x, y, t)dy = 1. 

We notice that, by the assumption, A(x) > 0 is a regular function on dVt and all 
other coefficients are constant. By definition u + is the solution of the system 

(d t - D+ A) u+ = 0, leOcK", 

U + \t=o = 1 ) 

du- 
\ dn 

which can be reformulated for v = 1 — u + 


^ Ian, A > 0, 


“+|an — M- 


(d t - D + A) v = 0, xeflC M n , 
'£>U=o = 0, 


(1 — h)|an 



D„du\ 
A dn ) 


an, 


where 0 < < 1 for all t. Moreover, as 0 < v < 1, it follows that 


0 < 



d U- A 
A dn ) 


an < 1 


and, as U- is increasing in time on dVt , then v is decreasing in time on dVt. Therefore, 
v < v , where v is the solution of the following problem: 


(d t — D + A) v — 0, x G C M n , 
v\t= 0 = 0, 
v|an = 1, 
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Thus, as in Ref. [29] (p.231 Lemma 12.7) for n = 2 , but now in R" (n > 2), we find that 
for the ball hi = U r (z) centered at z and of radius r , the solution satisfies as t +0 


v(z, t) < C 


/ \ n—Z 

ivm) exp 





with a constant C > 0 depending only on n (C can be explicitly obtained by the 
integration by parts in the generalized spherical coordinates in M n , where the coefficient 

/ \ 71—2 


I. 


+0 ° e '”' 2 ®"” 1 


corresponding to the leading term as t —» + 0 , appears from the integral 


l dw). Consequently (see Ref. [29] Corollary 12.8 p.232), for z £ int{Q} 


and t —y +0 we find 


v(z, t) <C 


dist(z, 


n —2 

exp 


/ dist( 2 :, 9h2) 2 \ 

V 4D + t ) ■ 


Then we immediately obtain Eq. (1171) by integration. 

For n = 2 we obtain directly the exponential decay in Eq. (fI71) for all e > 0. If 
n > 2, we still have the exponential decay for a small constant a > 0 depending only on 
e: 

0 ((7Ib)’’" e " 2/(4D+ ' , ) =0(e '“ / ‘ , ' 

Note that O ^e -e2 /^ AD + t ^ gives an exponentially small remaining term iff e = 2 \JD + t^~ 5 ° 

for a constant <5 0 > 0. For small enough <5 0 we have e = 0(y/4D + t ), also knowing that 
e > yJAD+t. 

So, for this e, we split hi in two parts: h2 e , the neighborhood of dVt such that 
dist(x, dVL) < e , and £l\£l e . For all F C hi \ , dist(F, dQ) > e > 2 y/D + t , we have 


/ (i - / G(x, y, t)d y'j dx = O [e c ( F )/ tS(F) 'j ? 

where c(F) and 5(F) are positive constants depending only on the distance between F 
and dVt and the dimension n. 

To complete the proof of the second statement, we first find that 


N(t) = / u_(x,t)dx 

J R n \n 


JR n \n Jn 
= Vol(ft) 


G(x,y,t)dydx = 



Js. n Jn 

G(x, y, t)dydx. 


G(x,y,t)dydx - / / G(x,y,t)dydx 
Jn Jn 


n Jn 
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For f2 = 0 e U (hi \ fi e ) we can write 



G(x,y,t)dydx = 


n Jn 



n € Jn Jn\n e Jn 


G(x, y, t)dydx 



= voi(o \ a) + 

in e Jn 

= Vol(fi) - Vol(a) + 


G(x, y, t)dydx + O [e t 7 
G(x, y, t)dydx 




'Q e Jn e 

G(x , y, t)dydx + O ( e - * 7 


Moreover, 



n e Jn\ n e 


G(x,y,t)dydx = 


'Q e Ja.\n e 
= Vol(fi\a)+ 0 (e - ^ ) - 



G(x, y, t)dydx — 


n Jn\n t 


'o\n e Jn\n e 


G(x, y, t)dydx 


' 0\f2 e J Q\Vt t 


G(x,y,t)dydx 


and since 


' o\n e J f!\f2 e 


< 


' f2\O e J Q 


G(x, y : t)dydx — Vol(f2 \ Q e ) 

G(x, y, t)dydx — Vol(f2 \ f2 e ) = O (e~i 7 'j , 


we conclude that 


and finally 



G(x,y,t)dydx = O [e t 7 


n e Jn\n t 


N(t ) = Vol(a 



G(x,y,t)dydx + O (e t 7 J , 


f2„ J fl e 


that completes the proof. □ 

A variant of Lemma |T] for n = 2 can be found in Ref. CD, where the heat localization 
near the boundary is also called by the principle of “not feeling the boundary”. In addition, 
we can consider the case of the distinct parts of the boundary: 


Corollary 1 Let X and Y be different closed parts of dLl such that dist(A",R) > 2e, 
where e = 0(\/t) > 2 y/Df t. Let U r (X ) = {x G R n | d(x,X) < r} be the open 
neighborhood of X of size r > 0 . Consider u + and u + as the respective solutions of the 
following systems: 


d t u + — D + Au + = 0, x G C M n , 
U+ U=o = 1 , 

D- du- 
\ dn 


^ \ash A > 0, 
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d t u + — D + Au + = 0, x G Cl C M n , 
u + \ t =o = 1 , 

„ I / D- du_ ^ ! 

u +\dvmu(x) = ( u ~ y Qff J I<9nnt7(x)> A > 0 

= 1 , 


u + \ 


dCL\(dflrU(X)) 


where U(X) is an open neighborhood of X of a radius strictly greater than 2e : C/ 2e (JC) C 
U(X). 

Then there exists 6 > 0 such that 


'Ue(X) 


\u + — u + \dx = O (e ** ) . 


Moreover, if u + is the solution of the system: 


d t u .|_ — D + Au + = 0, x G C M n , 


U + \t=o — 1, 

D_du\ 
A dn ) 

w+|an\Y = 1, 


u +\y — 


u 


Y, 


A > 0 


then 


'Ue(X) 


(1 — u + )dx = 


in\u e (Y) 


(1 — u + )dx = O [e 


The proof of Corollary [T] follows from the proof of the first statement of Lemma [TJ 

Note that the continuity of u on A (see Theorem ED and the localization of the heat 
content near the boundary allow one to consider mixed boundary conditions: 


Theorem 3 Let Q be a bounded (e, 8) -domain (see Section [||) with a closed connected 
d-set boundary dfl = Tq U T a U Too . Let A G C(To U T A ) such that 


f 0, x G T 0 , 

A(x) = < 0 < f(x) < oo, x G T A , 
I +oo, x G Tqo 


and e = 0(\/t ) > y/AD + t . We assume that the connection between different types of 
boundary is performed in the continuous way (see Theorem U|) such that the solution u 
remains continuous as a function of A . 

We split the e -interior Minkowski sausage of dfl into disjoint subsets 

a = u n r e x u 

such that each subset is contained in the e -interior Minkowski sausage of T (T C 
dTl). Then, for 8 > 0 from Lemma El the heat content of the problem (E]j-([^p, 

N(t) = j (1 — u(x, t))dx — j (1 — u(x, t))dx + 0(e~tz), 

Jfl Jfl e 
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can be found as a sum of two heat contents: 


N(t) = / (1 — u(x, t))dx + / (1 — u(x, t))dx + 0(e t 3 "). 

Jo r C Jo£°° 

In order to locally approximate the solution of the problem ([Tj)—(SJ) by considering the 
problem with coefficients frozen on a fixed boundary point, according to Corollary [T] we 
also obtain the following proposition: 

Proposition 2 Let a be a fixed point of the boundary <90 and let define 

= u u (a) n (O e U 0_ e ) for l G N, (19) 

where Ui e (a) C W 1 is a ball of radius le centered at o, e is defined in Corollary [B Let 
4> a G C°°(B 4ete (cr)) be a smooth cut-off function with a compact support on B 4e>e (o) : 

( 1 x G B^fa), 

4> a (x) = < a smooth function 0 < r] < 1 x G B 4e)e (a) \ B 3£:C (a), (20) 

{ 0 iGH \ B ie>e (a) 

If u is the solution of the problem (EP-@; then <f> a u is the solution of the following 
problem: 


d t - D ±A (< Ku ±) 


-(1 - u±)D±Acj ) (r ) 
0 


(4>crU±) U =0 = ln(x)<l>a(x), 
d((f a u-) 

D -—x-= \{x)<p a (x)(u- 


«+)bo, 


d{<j> a u +) d(<j> a u-) 

—x-ao — —x- 

on on 


x G B 4e}£ (a) \ B 3e>e (a), 
elsewhere in hi, 


( 21 ) 

( 22 ) 

(23) 

(24) 


Therefore, there exists 5 > 0 such that 


|u — (fa-uldx = O (e 


J B 2e ,eO 

and if (j) a u a is the solution of the problem / fUTD -pm) with frozen coefficients in the bound¬ 
ary point a, then 

(25) 


' 0\B e ^{a) 


r(1 — u a )dx = O [e t 7 


4 Local coordinates for a regular 8Q.eC 3 

In order to prove Eq. (J6]) for a large class of (e, 5) -compact connected domains hi in 
M n , we first prove it for the case of domains with regular boundary <90 G C°° or at 
least in C 3 . As hi is compact, for all types of connected <90, the volume of Id is finite 
and, therefore, the volume of the e-neighborhood of <90 in O is also finite and can be 
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approximated by a sequence of volumes of Minkowski sausages with regular boundaries 
(the same argument was used in Ref. [TT] p.378). 

Let us consider the regular boundary dfl G C 3 . 

Given a positive e > 0 provided in Lemma [U we denote by Q e and Q_ e the open 
e-neighborhoods of dVt in hi and in R n \ hi, respectively. 

According to Eq. (fT 8 l) and the regularity of the boundary dVt , we can decompose 
Q e U dVt U h2_ e = |_j[ =1 B i e (I is a hnite integer because fl + U is a compact domain) 
in such way that on each B i e it is possible to introduce the local coordinates. In addition, 
we assume that for all i — 1,..., / there exists u* G dfl D such that C i? 2 e,e(o'i) 
(see Eq. (fT9l) for the dehnition). Due to Proposition O the last assumption ensures that 



For all i we perform the change of the space variables (aq, ..., x n ) G B i t to the local 
coordinates (6i,, 0 n -i, s ) by the formula 



(26) 


where x(9i, ..., 0 n _i) G dfl and x , x and n are the vectors in M n such that 



is an orthonormal basis in M n . 

In what follows we denote Bi t D by Oj i+e and Bi, e D D_ e by ilj _ e respec¬ 
tively. In each of two regions, Qj j+£ and 0^_ e , the change of variables (xi,... ,x n ) h-)- 
( 0 i, ..., 0 n _i, s ) is a local C 1 -diffeomorphism. 

In local coordinates d£l is described by s = 0 . 

Thus, Eq. (fT 8 l) becomes 



(27) 


Denoting 6 = (0 l7 ..., 9 n _ i), the integration domain in (1271) becomes 


= {0 < s < e, 6 G dQ fl f\ +e }, 


which is actually a parallelepiped neighborhood (dQ fl Dj. +e )x] 0 , e[. 
For this change of variables we have 
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and therefore, using twice the integration by parts and the notations 


71—1 

iJM)i=n(i-^) ( 28 ) 

Z=1 

for the Jacobian and k t = ki(9i, ..., 0 n _i) of the principal curvatures for dfl curving 
away the outward normal n to 0 f 2 like in the case of the sphere, we find that for all test 
functions 0 = ( 0 + , 0 _) G V\ Bi>e 


'B, 


d t u | J(s, 0) | 0 dsd0i • • • d6 n -i 

^ \ n—1 

nn , \ 

D+\J(s, 9)\ 

D-\J(s,9)\ 


d_ 

ds 


— € 


d_ 

ds 


,du + 

ds 

, du _ 


E 


ds 


d / D + \ J(s, 9 )| du + 

Tl dfi V (1 - sk i) 2 "i 

d f _D_| J(s, 0)| du _ 


Z =1 
72 — 1 

E 


— as i V (i - f*i ) 2 as, 


0 + dsd0i • • • d9 n - 1 
c t>-dsd9i ■ ■ ■ d9 n -1 


+ / A(0)(w+ — m_)(0 + — 0_)<Z0 = 0. 


Is= o 


The regularity of the boundary ensures that the principal curvatures ki(6 ) are at least in 

c^dnndB^). 

Therefore, the problem (JT])—(J4j) locally becomes 


d 


d 2 


72—1 


d 2 


72 — 1 


m u+ ~ D+ \ te^ + Etw )“+ = fl +E — 


ski(9) 


ds 2 


72—1 


7=1 

72—1 


(90? 


1 + 


—' 1 — s/c,(0) y 1 — ski(9) J d9 2 


d 2 u , 


--D+ E s E TV 


k?(6) \ du + 


. Z— 1 


— J 1 — ski(9) ) ds 


i=l 


D a 


n—1 


Y — 

on. 


5 / | J(s, 0)| A 0W+ 


|J(s,0)|^00, V(l-^(0)) 2 ; 00* 


-, 0 < s < e, 0 G (<9f2 D 0* i+e ) 


(29) 


d 


d 2 


72—1 


d 2 


72 — 1 


( Oe2 / On^ 


0S 2 


'n—1 


*=1 

n—1 


de 2 


M_ = 


*-E- 


ski(9) 


1 + 


— 1 - ski(9) V 1 - ski(6 ) y (90 2 


0 2 u 




fc 2 (0) \ 0w_ 


, Z— 1 


—^ 1 — ski(9 ) y <9s 


L>_ 


72—1 


E on. 


<9 / | J(s, 0)| \ (9u_ 


17(5,0)1^,90, V(i-^(0)) a ; 50, 

^+|t=0 1; Bj— |t=0 0, 




0W_ 

ds 


s=— o A(0)(w_ ) | s=0, 


ZA 


0 -u 


+ 1 


ds 


|s=+o — D. 


du _ 
0 s 


I s=—0- 


-, -e<s<0, 0G(0O nh, i+£ ), (30) 

(31) 

(32) 

(33) 


19 





























We emphasize that the problem (j29|) - fl33|) should be considered as the trace of Eqs. ©-o 
on B i t in the sense of the problem (jZD)~(123D with (f)g i = 1 on B ie . 

Therefore, we can rewrite (12711 in new coordinates and use the parallelepiped property 
of Qj )+e in the space of variables (s, 9) : 

1 r 

N(t) = / (1 — u(s, 9, £))| J(s, 9)\dsd0 + 0(e~t s ) 

l=\ 

= ^ j d6 f ds(l — u(s, 9, £))| J(s, 6)\ + 0{e~t s ). 

i = \ JdQiC\Q.i ,-|_g J[ 0,e] 


Since this local representation holds for all i (the form of the problem (l29]) - ([33]) is the 
same for all i ) and Xli=i Ion no + 99 = J gQ d6 , we can formally write 

N(t) — f d6 f ds(l — u(s,9,t))\J(s : 6)\ + (34) 

Jdn J[Q,t] 

where u is the solution of (l29 |) - ([33j) in ] — e, e[x<9fl in the local sense, as explained 
previously. 


5 Approximation of the heat content by solutions of 
one dimensional problems (for a regular boundary) 


We denote by G(si, 9\, s 2 , @ 2 , t) the Green function of the problem (l2f)l) - ([33j) in <9fix] — 
e, e[. Let us fix a boundary point (0, 9 0 ). 

We denote by G e ° the Green function corresponding to the following constant coeffi¬ 
cient problem, considered as a local trace problem, i.e. in the sense of the problem (1211) - 
(J23J) with (j) e . El on B ie : 


Next, let 


d 


d 2 


n—1 


d 2 




dt 

d 


ds 2 
d 2 


i=l 

n—1 


se'f 
a 2 


< s < e 


m u - ~ D - 1 ^ + E* 1 = °. - £ < s < 0 


ds 2 


i—l 


DO 2 


u + \ t= o — 1 , U-\ t= o — 0 , 

dzi— 

D —— I s=—o A(6 l 0 )(n_ , U-)_)| s= o, 


D a 


ds 
du+ 


ds 


|s=+o — D-- 


du _ 
ds 


s=-0- 


(35) 

(36) 

(37) 


(38) 


^R"( S 1, ^1) S 2, 92 , t) 


l{si>o}G ? ++( s i) 9i,S 2, 9 2 , t ) + 


l{si<o}G ,0 °_ ) _(si, 6i, s 2 , 0 2 , t ) 
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be the Green function of the constant coefficient problem in the half space, explicitly 
obtained in [B] Then, according to Ref. |30| p.48-49, due to Varadhan’s bound property 
of Green functions, in U e { 0, 0o) the difference between the Green function 4>g 0 G 9 ° of the 
problem (l35li (1381) and the analogous Green function in R n , , is exponentially small: 

\{J>e 0 G do - |c/ e ( 0 , 6 >o)x?y e ( 0 , 6 »o) I = O (e ^ 


Therefore, following the ideas of McKean and Singer[30] (p.49), we approximate G by 
the Green function G e ° with the frozen coefficients on (0, 6 0 ) , whose replacement by G^ n 
yields only an exponentially small error. 

For an abstract operator Cauchy problem 

8 

—u — Au = TZu, (39) 

u\t=o = Uo 

the solution u can be found by the Duhamel formula 

u(t) = e~ tA u 0 + [ e- (t - T) ^u(r)dr. (40) 

Jo 

Therefore, by the Duhamel formula, locally, we have the following infinite expansion 
u+(s,0 o ,t)= I d0 1 ds 1 G^ + (s,0 o ,s 1 ,0 1 ,t) 


+ dr d6idsiG e ° + (s, 9 0 , s\, 9i, t — r) • 

Jo Jn e 

■ f dd 2 dS 2 'R-G+ + (si, 9i, S 2 , & 2 , t) 

Jn e 

+ f dr f d9\ds\G e ° + { s i 0 O , s\, 0 3 , t — r) ■ 

Jo Jn e 

■ f dri [ d9 2 ds 2 TlG%(s 1 ,9 1 ,S2,92,T - Tx) ■ 

Jo J?L t 

• f d0 3 ds 3 KG%(s 2 ,02 , S 3 , 03 , Ti) + .. . + O (e~* 


where the operator 1Z is dehned by 

1Z = TZ Sl (si, 9\) + TZg 1 (si, 0i), 
8 


/ n—1 


n —1 


TZ s (s , 6 ) = R(s, 0)^ = -D + [J2 h{9) + s 


^{0) 


8 


, i=1 


^ 1 — skAO )/ < 9 s’ 
1=1 x 7 / 


(41) 

(42) 

(43) 


n —1 


K„( S ,0) = ZH m 


1 + 


^ 1 - ski{0) V 1 - ski(9) J 891 


8 2 


D 4 


n—1 

£ 


8 


|J( S ,0)| \ 8 


\J(s,9)\^89 i V(1 -sHO)) 2 ) d9i ' 

We substitute Eq. (j4Th into Eq. (j34l) with 9 = 9q and prove the following theorem: 


(44) 
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Theorem 4 Let 


u .|_, 0 < s < e 
U ^ m_, — 6 < s < 0 

be the solution of the one-dimensional problem 
d d 2 

—u - D±——u = 7 Z s (s, 9 0 )u - e < s <e, 9 = 9 0 , 

at os 2 

^|t=0 lo<s<e(s); 

du 

D_ — \ s= _ 0 = \(9 0 )(u_-u + )\ s=0 , 


D 


du A 


~\~ r\ <S —+0 

as 


—4-n — D- 


du- 

ds 


|s=—0) 


(45) 

(46) 

(47) 


obtained from l2£h) - l3^) setting 9 = 9 o (TZ s (s,9o) is given by (f IjBJ). Then the heat 
content N(t) , defined in a, satisfies 


N(t) — f d9 0 [ ds(l-u(s,9 0 ,t))\J(s,9 0 )\ = I 0 < A < oo 


Jon j[o,£] 

If all principal curvatures of dLl are constant, then 


0(t 2 ), A = oo 


(48) 


N(t)= / dO o / da (1-u(a,0o>O)k(Mo)l+0(e <*)■ 

JdQ. J[ 0,e] 

Moreover, if u hom is the solution of the homogeneous constant coefficients problem 


d 2 

d t u - D±—= u = 0, -e < s < e, 9 = 9 0 , 
os z 


(49) 

& 1 1=0 lo<s<e(^)? 



du 

D_ | s= _o = A(0 o )(tt- tt+)| s =o, 


(50) 

du+ _ du_ 

+ ds ls=+0 " |s= -°’ 


(51) 

f d0 a ! da(l-u fc “(a,#„,i))U(a,#o)l = ( 

Jdfl J[0 ,e] l (J \ t h 

0 < A < oo 

A = oo 

(52) 


then 

m 

From E we get 

G e °(s l ,9i,S2,9 2 ,t) — 1{si>o,s 2 >o}G ! + 0 +('Si, 9i, S‘2 , 9 2 , t ) 

+ l{si<o,s 2 >o}G ,ft) + (si, 9i, s 2 ,9 2 , t). 

Due to Eq. (1341) . we need to know only G°fi + 

Gq? + (si, $i, s 2 ,9 2l t) = (/4°(si, s 2 , t ) — f+(si, s 2 , i)) A (fix, 9 2 , D + t), 
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with notations 


h+(si, s 2 , t) = 


exp 


+a(A, 0, 0 O ) exp 


a/47t D + t 

(si + S2) 2 


(si — £2)^ 

4 D:t 


4 D+t 


fl°(s 1 ,s 2 ,t) = b(X, 0,0 o )^j^- exp + s 2 ) + A (9 0 ) 2 a 2 t 


D + ^WDT 




where 


= , A (0 O ) = 00 , 


1 , A(0 O ) < 00 , 

a{\, 0, 9 0 ) = ^ y/pf-y/pZ 
^/pZ+^/pZ 

1, A(0 O ) < 00 , 

0 , A(6 I 0 ) = 00 , 

and K(9i, 9 2 , D±t) is the heat kernel in M n_1 : 

1 


&(M,0o) = 


K(9 1 ,9 2 ,D±t) = 


(4nD±t) 2 


— exp 


\9i-9 2 \ 2 

4 D±t 


Since 


N(t) = / (1 — u e (s,9,t))\J(s,9)\dsdd + 0(e ^) 


= Vol(a 



G(s, 9, si, 0i, £) | J(s, 0)|dsd0dsid0i + 0(e ** 


n £ </n £ 


(53) 


(54) 


(55) 


in what follows we use P(t) for the notation of the principal part of N(t): 

P(t) = I f G(s,9,si,9i,t)\J(s,9)\dsd9dsid9i. 

Jn t </n £ 

To prove Theorem [4] we need the following Lemma: 

Lemma 2 The principal part P(t ) of the heat content for the solution of the system 
Q5j) . defined in Eq. K51 1) . is given by 


P(t) — / / G e ? + (s, 0Q; s i ; 0i, t)|J(s, 8 0 )\dsd9 0 d9idsi + G^AifP 6 ^ 


(56) 


(57) 


n £ Jn e 


+n e e °)G d z + + G e zj(n 9 s ° + n e e 0 )G 9 z4{n 9 s o + n 9o )u+ + o( e ~*), (58) 

with notation 

G 9 °J(n 9o + n 9 °)G e ° + = f dT [ dsd6 0 \J( Sl 6 0 )\ [ G+ + (s, 0 O , s l5 0 1; t — r) • 


'o ./n £ 


(ft£(s l5 01) + ftJ°( Sl , 0i))G® o + ( Sl , 01, s 2 ,02, T)d0 2 ds 2 d0idsi. 


(59) 
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Moreover, the following identities hold 


M(t) — 1 f G+ + (s,9 0 ,si,9i,t)\J(s,9 0 )\dsd9 0 d9 1 dsi 

Jn e Jn e 


= ) 

f d9 0 [ dsids (h 0 _f!(s,si,t) - fl°(s,s 1 ,t))\J(s,9 0 )\, 

an d[o,e} 2 

(60) 

G 6 ° + pi e °G e ° + = G\^R e e°G^4n e e °G\ = ... = 0, 

(61) 

Jon 

d9 1 K(9 0 ,9 1 ,D + t) 


= J 

( d9 1 [ d9 2 K(9 0 ,9 1 ,D + (t-T))K{9 1 ,9 2 ,D + r) 
an Jan 


= . 

• • = 

(62) 


Proof. Formula (1551) is the direct corollary of the Duhamel formula (see (I4U1) and (TUT) ). 
Let us start to prove (l60lh 
Indeed, we find that 


M(t) 



f+(s, si, t))K(6i, d 2l D + (9 0 )t)\J(s , 9 0 )\dsidsd9id9 0 




ddoddi .„-i exp 

(AllD+t) 2 


( 16 > 0 - 6 >i| 2 \ 

V 4 D+t J 


lan(9o)^-dn(9i)$(9o, t), 


where 


®(9 0 ,t)= / dsds 1 (h e f(s,si,t)-fl°(s,si,t))\J(s J 9 0 )\. 
J[ 0,e] 2 

With the change of variables 9 1 H> v = , M(t) becomes 


(63) 


M(t ) = 


f e-l^l 2 

/ n — 1 

i^- 1 ir- 1 7T 2 


lan( 6 | o)l 9 n+ yj 5 ^ 7 ^( 6 l o) 4 ) ( 6 'o, t)dvd9 0 . 


By our construction, 


9 0 G dQ and 9 X = 9 0 — \JdD + tv G <9f2, 


that implies 

lao(^o) - 1 ^( 0 o)l af2+x /^ w (^o) = 0. 

It can be interpreted in the following way: if we take a point on the boundary and move 
it along the boundary, we obtain another point which is still a boundary point. 
Consequently, we find (1501) 


M(t) 



^——rtdQ.(9o)$(9o, t)dvd9 0 = 

7T 2 



H9 0 ,t)d9 0 , 
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which also implies the first part of (j62j) : 


d9 1 K(9 0 ,9 1 ,D + t) 


'an 


lan( 0 o)- 


Let us now prove that in the computation of P(t) all terms containing the derivatives 
over the transversal variable 0 vanish. 

For all terms in (091) containing a derivative over 9 \, we calculate (see (|Hj)) 


n ei K(9 1 ,e 2 ,D+t) 

D + siki(6i 


n— 1 


1 - sifcj( 6 >i) V 1 — siki{9i) 


( 0 j - eif 


2 D + t \ 2L> + t 


-l)K(9 1 ,9 2 ,D + t) 


n —1 




d f D + \j( Sl ,9 1 )\\(e\-9i) 


\J(s 1 ,9 1 )\^d9\ V(1 -sMOi)Y 
Let us prove Eq. (RTTlh noting that 


2 D+t 


-K(9 1 ,9 2 ,D + t). 


G e °4n dl G e ° + = / dr / dsd9 0 \J(s,9 0 )\ / G 0 + ° + (s, 0 O , s u 9 u t - r) • 

Jo Jn t Jn t 

■ / n ei G^ + ( Sl , 0 i, s 2 , 0 2 , r)d 0 2 ds 2 dM Sl . 
d Q e 

We can schematically rewrite G e + + $'R.Q 1 G e + + in the following form: 

GjJRd.Gjv = 

= f dr f de„ f Mif d8 2 K(Oo 1 0 1 .D + (t — t))1Zo 1 K(9 1 ,92, D + t) 
Jo Jr™- 1 J R"- 1 Jr™- 1 

dsdsids 2 0(s, si, s 2 , t, t, 9 0 ). 


■ 1 ao (' 0o) 1 an ('J i ) 1 an (' 02 


'[ 0 , e ]' 


With the change of variables involving 9q : 

60 -el 


9 1 = 


^y/ D +( t ~ T ) 


, 0 1 = 0 o - 2 v / J D + (f-r) 0 ~ 1 , 


n 1 _ /q2 

02 = 0 2 = 01 - 2 Jd^9 2 , and 

2 /D 7 W V + ’ 

02 = 0q — 2y/ D + (t — r)0i — 2^/ D + t9 2 , 


so 


(64) 


( 65 ) 


and since for all 0o G <912 


W0o) - lan(0o)l an+2 ^ £ , +(t _ T) gi 1 (0o) ' 
'^an+ 2 A /D + (t-r)ei+ 2 x /D + r 0 2 ~ ^ 
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we obtain the separation of variables on 9 2 from (9 0l Si,9i): 


G++tl^0iG'+ o + = dr d9 0 / dsdsids 2 0(s, s u s 2 , t, r, 9 0 ) 


'++ 

n—1 

n 

i=l 


d9\ 


,e-(«i) 2 


d.9 


an 

. e -d>\? 


'[ 0 , 6 ]: 


7r 


7r 


C\(2{9\Y-l)-C\9\ 


where C\ and C 2 are the functions of Si,9q,9i , but not of 9 2 , and consequently 

G^JUeG 9 ^ = 0 . 

By the same reason we have Eq. (16TT) . Changing variables 9 t to 9^ from (1641) (l65ll . we 
also obtain the last part of (1521) . □ Let us know prove Theorem U 

Proof. To End Eq. (jj5J) , we study Eq. (1551) using proved relations (jBDl) - (1521) . For 
instance, we have 


G e °4(n d s ° + n e g °)G\ = G e ° + pi e s °G d ° + 



As G\(s : 9 0 ,si,9i,t) = (h e °(s,s u t) - f+(s, si, t))K(9 0 , 6h, D + (9 0 )t) , we have 

G^JTZ e s °G% = (h°° - f^)Kpi d s 0 {h^ - 

Now we perform the change of variables (I64ll - (l65lh Since locally ki{9) G C 1 , then for 
all 9 0 e , for t —> +0 we can develop 

ki(9 0 - 2a/ D + (t — r)0\) = h{9 0 ) - \/k i (9 0 )2^D + (t-r)9 1 + 0(t - r). 


Consequently, by definition of 1Z sl (si, 9 1) in fl4,‘j|) . which is a composition of the operator 
of the first derivative by Si and of a multiplication by a function of the class C 1 on 9\ 
(locally, in the sense of local variables), we also have for t —> +0 


ft Sl (si,0i) = Tl Sl (si, 9)[1 + 0(t - r)] - V e Tl Sl (si,0)2y/D + (t - r)0i. 
As 2\jD + (t — r) J R d9je^^ 2 9j = 0 , we obtain 

G e °jn e s °G 9 ° + = f dr [ d9 0 [ ds\J{s,9 0 )\ / d Sl [ ds 2 


'an 


'[ 0 , 6 ] 


'[ 0 , 6 ] 


'[ 0,6 


(h+ — /+°)(s, Si,t ~ t) 1Z s1 (si, 9 0 ) [1 + 0(t — r)] ( h e ° - /J°)(si, s 2 , r) 

n—1 


n 

i= 1 


d9\ 




e v > 


d9\ 


-(C 2 ) 2 


e 


7T 


7r 


= dr d9 0 ds|J(s,6> 0 )| / ds x 


ian 


'[o,. 


'[o,. 


'[ 0 , 6 ] 


d s 2 (h e ° - fl°)(s,si,t - r)n si (s ll 9 0 ) [1 + 0{t - r)] (h 6 ? - /+°)(si, s 2 ,r), 
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from which it follows 


P(t) = / d9 o / d s dsi(h e ^(s,si,t) - /J°(s, s 1 ,t))\J(s,9 0 )\ 

Jdn J [ 0 ,e] ./[ 0 ,e] 

+ [l + 0(t)] f dr [ d0 o ( ds|J(s,0 o )| [ ds^h 9 ? - f+)(s, si, t - r) 


/an a[0,e] 


f [o, e ] 


'[o,d 


ds 2 n si ( Sl ,9 0 )(h d + ° - f e + °)( Sl ,s 2 , 


T 


[1 + 0(t)} 2 / d9 0 (h e ^ - f d + °)Wl s {s,9v)(h e ° - f e + °)$1l s {s,9 0 )u A e 


( 66 ) 


I an 


We notice that the solution u(s,9o,t) of the one-dimensional system (H5]) - (I47|) is given 
by 

u(s,9 0 ,t)= [ ds 1 (h e £(s,s 1 ,t) - /J°(s,si,t)) 

J[0,e\ 

+ f dr [ dsi(/4° - /J°)(s,Si,t - r) f ds 2 K 81 (s u 9 0 )(h e _£ - /+°)(si, s 2 , r) 


/0 J [0,e] 


'[o,d 


+ (/ 4 ° - 

To obtain (j48l) of Theorem |4] from formula (I66[) . we estimate 


NN 2 (t) = 0{t) I dr / d9 0 / ds\J(s,9 0 

Jan J[ o,e] 


'[0,6 


0 -JO\L o |u,e| 

dsi(h+ — /+°)(s, s\,t — t) I ds 2 ^ si (si,0 o )(/i+ - /+)(si,s 2 ,r). 

•/[o,e] 


In fact, from f|52p . proven in what follows, it holds (see (1551) for the definition of iViV 1 (t) 


NN 2 (t) = 0(t)NN\t ) = O(t) 


O(tz), 0 < A < oo 0(t 2 ), 0 < A < oo 


O(t), A = oo 


0(t 2 ), A = oo 


To conclude, we note that if all principal curvatures kj{9) on are constant, then 
for all 9 e dVt 

Tl s {s,9) = 77 s (s,6» 0 ), 

and thus 

N(t) — j d9 0 f ds (1 — u(s, 9o, t))| J(s, 0 O )| + 0(e~t*). 

Jan J [o,e] 

To show (152|) . we need to estimate 


NN‘(t) = £ rJvttft.rf+t • • • (i - fold), 


( 68 ) 


Z=1 


where 


r e + ° + = (^° - /+) 
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More precisely we want to prove that for all j > 1 


1-i-j 

t~y/AD+t), 0 < A < oo 


\NN j (t)\ < C 


y/AD+t), A = oo 

Due to Lemma [21 we start with (see 


(69) 


NN 1 {t)= dr d9 0 / ds| J(s, 0 O )| 

Jo Jan J[ o,e] 

dsi(h+° - f*>)(s,si,t - r) [ ds 2 TZ sl (s 1 ,d 0 )(h e ° - fl°)(s 1 ,s 2 ,r). 


J[ o.d 

Therefore, we have to estimate four terms: 


'[0,6] 


where 


NN l (t) = ^ 

3 =1 


MM\(t)— dr (W 0 / ds\J(s,9 0 ) 

Jo Jan J[ o,e] 


'[ 0 ,' 


dsih+°(s,si,t - r) / ds 2 ^ Sl (si, 6 l o)h^ 0 (si,S 2 ,'r), 


'[0,6 


MM 2 (t) — — dr d9 0 / ds|J(s,# 0 ) 

Jo Jan J [o,e] 


'[0,6] 


dsi/+(s,si,t - t) / ds 2 n si (si 1 9 0 )h e °(s 1 ,s 2l T), 


'[ 0 , 6 ] 


MM 3 {t) — — dr d9 0 / ds|J(s,# 0 ) 

Jo Jan J [o,e] 


'[o,. 


dsih e °{s,si,t - r) / ds 2 7^ Sl (si,6>o)/+ 0 (si,s 2 ,r), 


'[ 0,6 


MM 4 (t) = dr cW 0 / ds|J(s, 0 O )| 
Jo Jan J[ o,e] 


ds 1 fl°(s,s 1 ,t - r) I ds 2 'Jl sl (s u 9o)fl°(s u s 2 ,r) 


'[ 0 , 6 ] 


We aim to approximate TZ Sl (si,9 0 ) = R(si,9 0 )d Sl from Eq. (1431) near the point ( s,9 0 ) 
For t —y +0 and 0 < Si < e = 0(y/i ), we End that 


1 - Si/q(6> 0 ) 


— 1 + siki(9o) + O(si), 


which gives 


' n— 1 


n— 1 


R(si, 9 0 ) — —D + j ki(0o) + si k'i(9o) + 0(&l) 


, 2 — 1 


i=l 
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Introducing the notations 


(si =b s 2 ) 2 
4 D + t 


c = 


1 


C± = si ± s 2 , 


I Si ±s 2 {t ) = exp 


we find 


8nDl(0,9 0 )y/(t - t)tz ’ 


h+°(s,si,t - t)71 Si (si, 6 0 )h 9 ^(si, s 2 , t) = —cR(si, 9 0 ) (C_[/ s _ Sl (i - r)/ Sl _ S2 (r) 
+a(A, 0, 9 0 )I S+S1 (t - t)/ s1 _ S2 (r)] 

+a(A, 0, 9 0 )C + [/ S _ S1 (t - t)I Si+ 82 (t) + a(X,0,9 Q )I 8+8l (t - r)/ Sl+S2 (r)]). 

We now change Si to z\ and s 2 to z 2 by the following change of variables: 

• for P SlTS2 (r) : z 2 = ^dIt and s 2 = ±Si =F 2 ^/D + tz 2 , 

• for P sTsl (t -t): Zi = 9 ^/^_ r) and si = ±s =F 2 y/D + (t - t)z x . 

Let us notice that t is a constant parameter and, as r takes its values between 0 and t , 
hence, z\ and z 2 are in M + or M. But at the same time 2 \jD + (t — r)z\ — s ± Si and 
2y/D + rz 2 = Si ± s 2 are bounded to the interval [—e, 2e] and hence are of the order of 
O(yft') . In what follows, we suppose that r and t — t have the same order of smallness 
as t: 

0(t) = 0 (t) = 0(t — t). 

Therefore, for 0 < si = ±s =F 2 1 / D + (t — t)z\ < e we have 

P sl (s u 6o) = [4>(9 0 ) t^(s,z 1 ,9 q )]-^-, 

where 

n— 1 

<j>{9 o ) = -D + Y,ki(0o), 

i— 1 

n —1 

■0(s, z 1} 9 0 ) = (s - 2a/ D + (t — t)z 1 )D + Y k K°o) + 0{t). 

i =1 

If we develop R in the neighborhood of (0, 9q) , we find 

R(s±, 9 q) = R( 0, 9q) + 0{Vt) = <i>(9 0 ) + 0(Vt). (70) 

For A = 00 on dfl , we simply have 

MM 2 (t) = MM 3 {t) = MM A {t) = 0, 


and 


\NN\t)\ 



= \MMi(t)\ = \h e ^n s h e _°\ 

dr ~^= [ d0 o f ds|J(s,0 o )| 
V T Jdn Jo 


< Cy/t/i(dQ, \JAlD + t). 
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By iteration of the proof, we show that liViV- 7 (t)| < Ct^n(dfl, y/AD + t) for j > 1 
Now, for A < oo , 


MMi(t) = h 9 nn s h e ° = - / dr 


dOr] 


7T 


dz 2 (j){0 o )z2e 22 


^ Jon VD+ 

y2xi(s,z 1 ,z 2 ) - [ d T^-= [ 
JO V T Jd : 


ds|J(s, 0 O )| / d^e Zl 


7 r 




ds\J(s,e 0 )\ 


d^ie" 


2 

dz 2 ijj(s, z\, do)z 2 e~ Z2 




„ 2 — 1 


2=3 


Here for v + (zi) = 2 y/ D+(t - r)^il R +(^i), v(z 1 ) = 2^D + (t - t)z 1 , w + (z 2 ) = 2v^t£ 2 1r+(z 2 ) 
and u>(z 2 ) = 2 ^D + tz 2 


Xi(s, z u z 2 ) = l [0 , e ](s)l[o, £ ](s - v( 2 i))l[o, e ](s - v(z x ) - w(z 2 )), 

X2 (s, Z lf Z 2 ) = l[0, e ](s)l[0 i£ ](s - u(«l))l[0,e](-S + v(z{) + W + (z 2 )) , 
Ah(s, Zl, Z 2 ) = l[0, e ](s)l[0, E ](-S + V+(^i))l[0, E ](-S + V+{z{) - w(z 2 )), 
X 4 (s, Z U z 2 ) = l[0, E ](s)l[0, E ](-s + w + (z 1 ))l[ 0i£ ](s - v+(z-i) + w + (z 2 )). 

Considering two formulas: 


0 < V < £ 

0 < s < v 

(71) 

— £ < V < 0 

£ + V < S < £ 

0 < V < £ 

0 < s < v 

(72) 

£ < v < 2e 

V — £ < S < £ 


V = l[o, e ](s) - l[o, e ](s - v)l[ 0)E ](s), and C = l[o, £ ](s)l[o, E ](-s + v+), 
we find that 

V ± o <*=» 

Ct^O « 

ft means that for 0 < v = v + < e, it holds 

V(s) = C(s) = l[o,«+]( s ) 

and for — e < v < 0 and e < v + = v + 2e < 2e it holds 

V( S ) C( 5 ) l[u+— e,e]( s ) l[e+ti,e](' s )- 

Consequently, we found the formula 


l[ 0 , e ](s) - 1 [0,e]( s - v)l[ 0 , E ](s) = l[ 0 , E ](s)l[ 0 , E ](-S + V+), (73) 

from which it follows 


Xl = l[ 0 , e ](s)l[ 0 , e ](s -v) - X 2 , X4 = 1[0,6](s)1[0, e ](-s + v + ) - X3- 
Therefore, we have 

4 

^2 Xi(s , z 1 , z 2 ) = l[ 0 , £ ](s) and Xi + X 2 ~ X3 ~ X4 = l[o, £ ](s) - 2 • l [0 , w+ ( 21 )] (s), 
2=1 
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which are independent of z 2 • Since 


z 2 e Z 2 dz 2 = 0 and / (2 z 2 — l)e Z 2 dz 2 = 0, 


we obtain exactly 

\MM\(t)\ = \hpi s h\ = 0 . 
For MM 2 we find in completely analogous way 

A(0„) 


2 f* . A Jt — r 


MM 2 (t) = —= / dr-i UlUQ - 

V 77 Jo v r iao 


c/0n 


ds|J(s,0o)l 


dZ i e 2 A(0 o) az i A / t=?+A(0°) 2 a 2 (t-r ) + A(0 o )a A /t^7) 


d^2 


0(6 l o)^2e * 2 (x 3 (s, z i, 2 2 ) + X4(s, zi, 2 2 )) 

f V* - t f Jn A(0o) ' 7 


dr- 


7T 


dd()—— 
T Jan D.\ 


ds|J(s,0 o )| / 


/R+ 


^e 0 ) azi Vi=^+\(e 0 ) 2 a 2 (t-r) Erfc ( Zi + A (0 o )aVt^7) 


dz 2 


2 

-0(s, 21, 9 0 )z 2 e ~ Z2 (y 3 (s, 0i, 2 2 ) - x 4 (s, z h z 2 )). 


Since 


Xa(s,Zl,Z 2 ) + X4(s,Zl,* 2 ) = l[0,e](s)l[0, e ](-s + V+(*l)) = 1[0,«+(*!)] (s), 
X 3 (s, Z 1, 22 ) - X 4 (s, z 1, 22 ) = 2 y 3 (s, 2i, 2 2 ) - 1 [0,e] 0)1 [ 0 , e ] (~S + U+(2i)), 


the parts of MM 2 , which contain the integration over s on [ 0 , 2 y/D + (t — t)zJ\ , are 
equal to zero. In addition, for £ = 0,1 


d2i Erfc (^ + a {d 0 )aVt^T) 


< c. 


As %f> is of the order 0 (\/i ) and linear on 2i, and e = 0 (\/i) , we directly obtain 


\MM 2 (t)\ = 


2 ^ dr Xf-IfL f f &s\J(s, 6 a )\ [ dz, 

n 1 0 v r </ao *-'+ Jr J r+ 

2 VU 


. e 2A(e 0 )a~n^+A(@o) 2 « 2 h-r) Erfc ( Zl + A(0 o )aOf^T) 
• / d2 2 ^(s, 21, 6 l 0 )2 2 e _Z2 2x 3 (s, 21, z 2 ) 


<C / drA/r / d 6 0 / ds|J(s, 0 O )| < Ct 2 n(d£l, ^/dD + t). 
Jo Jdn Jo 


Since /i(< 9 Q, ^/ 4 £) + t) = Ca/^ for a regular boundary, then |MM 2 (f)| < Ct 2 . 
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To estimate MM 3 we find 


OsJ^s^r) = / -^ff + °( Sl , S2 ,T ) - 

v u + 

In our notations, using (T70li . we have 


^(0o) a £ e 0 /„ „ ^ 1 _ ( ( s i + s 2 ) 2 


D + yjnD + T 


exp 


4 D+t 


h e P(s, s u t - t)K > 1 (s 1 ,»o)/?(si,S2,t) = P ‘ *> + - P,+s ‘ 


\/4ir£>+(( - r) 


•(0(«o) + O(Vi)) ( ^=-f?°(si, S2 ,t) - A * 9o) * 


D + j¥D^ P ‘' +n 

Changing variables Si to Zi and s 2 to Z 2 , we obtain y 2 ± X 4 f° r t-he area of s , which 
gives intervals (linearly) depending on the values of z\ and Z 2 ■ Thus, we majorate s by 
e and estimate MM 3 : 

\MM 3 (t)\ <C\ f d Ty/r f d9 0 I ds|J(s,6> 0 )| 

Jo Jon Jo 


where 


d^ie z * / dz 2 (0(6> o ) + 0(Vt))f(z 2: T)\, 


f{z 2 ,r) = —^ 7 ==- exp (2X(d 0 )az2Vr + X(6 0 ) 2 a 2 T) ■ 


VD^ 

Erfc(z 2 + \{9 Q )ay/r) 


1 -z 2 

-.e 2 . 


We see that 


yJnD + r 
dz 2 {(j)(9 0 ) + 0{Vt))f(z 2 , r) 


< C. 


Therefore, we have 

|MM 3 (f)| < Ctn(dQ, y/AD+t). 

In the same way, since X 4 depends on z\ and Z 2 at the same time, we have 

dr \Jr(t - r) [ d9 0 [ ds|J(s,0 o )| 


\MM 4 (t)\ < C 


>dn 


p+oo 


dzie 2 X(e 0 )az 1 Vt=^+\(e 0 ) 2 a\t-r) Erfc (^ + X(9 0 )aVt^T) 


p+oo 


dz 2 ((f){9 0 ) + 0(y/t))f(z 2 , 


T 


< \JAD + t). 


By iteration of the proof, we show for j > 1 that 

\NN j (t)\ < Ct^^dQ, yJdD + t). 

□ 
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6 Relation of the heat content expansion with the in¬ 
terior Minkowski sausage 

Let us start with a heat problem with just a discontinuous initial condition. 


6.1 Particular case D + = D = const 

Lemma 3 Let 12 C M n be a compact connected bounded domain with a connected bound¬ 
ary <912 of the Hausdorff dimension d and u is the solution of the following problem: 


d t u — DAu = 0 x G M n , t > 0, 

(74) 

a 

o 

II 

-w 

”s~ 

(75) 

Then for t —>• +0 we have 


[•2 —z 2 

N(t)= j —^p(dQ,2\/lDtz)dz + o^t^~^j. 

(76) 

Moreover, it can be approximated by 


N(t) = / 3 n - d 2a /Dt) + o (t^'j , 

(77) 

with the prefactor 

[A - [ ZX£ * dz = 1 7 f ^ + 1 ,4) 

io A 20F ' V 2 ) 

(78) 


is expressed through the incomplete Gamma function. 
Proof. Let us prove formula (1761) . By definition 


N(t)= [ _ f G(x,y,t)t a dxdy, 

jR n \n JR n 

where this time G is the heat kernel in W 1 

G(x,y,t ) = (4Z97rf) _ 2 exp 


f_\x-yT\ 

V 4 Dt )■ 


Therefore, we have 

N(t) = Vol(12) 


./HJu, 


(dTrDt) 


\ \x — y\^ 

ne In ^tn(x)tn(y)dxdy 


Vol(12) — f —e ( f + 2\ /r Dtv)dx ] du 

J R™ 7T 2 \J R™ / 


_ P -H 2 

n ° 

7T 2 


(l n (x) - t n _ 2VWtv (x)) dx 


du, 


where 1 Q _ 2 y/Dtv( x ) = + 2 y/Utv) and the notation 12 — 2\/~Dtv means that 12 is 

shifted by the vector —2 \[T)iv e M" . 
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Let us firstly suppose that dfl is regular, i.e of the class C 3 . We see that for all points 
for which d(x, dQ) > 2V r Dt|M|, it holds (x + 2\fDtv) G hi. Thus, it follows that 
for e = 2 \/I}i||u||, 


In (a) (ln(®) - l n -2 VDiv( x )) = 0 for all x G hi \ 


Therefore, only x belonging to hl e with ||w|| < contribute to N(t) and we can 

write: 


N(t)= / — e 
Jm™ 7T 2 


-hi 2 


(ln £ (x) - t^_ 2VWtv (x)) &x 


d v + O (e ) , 


where the exponentially small error with a S > 0 is defined by the integral 

1 


—we ^ du. 


M > 


—£ _ 7T 2 

2 s/m 


Since dfl is regular, we introduce (see Section [JJ the local coordinates x = (9, s ) 
and thus have x(9) G clhl and x G hl e iff 0 < s < e. In this case, X 2 ^Dtv( x ) = 
lf 2 e (x) — 1 q £ _ 2 s/Dtv( x ) 7 ^ 0 iff x G hl e and x{9) — sn{9 ) + 2 \J~Dtv £ hi. Moreover, with 
the notation (v,n) for the Euclidean inner product of two vectors in R n , 


(. x{9 ) — sn{6) + 2 a J~Dtv) ■ n(6) = — s + 2 a /Dt(v, n). 


We deduce that 

X 2 y/m,v( x ) 7 ^ 0 iff s - 2 y/Dt(v, n) < 0 . 

Consequently, if (u, n) < 0 , as s > 0 , it is not possible to have s — 2 y/Dt(v, n) < 0 . In 
turn, if 0 < (v,n) then s g]0,2a /Dt(v, n)[. Considering only (v,n) > 0, we can define 
e = 2^Dt(v, n) and, since v = and x,y G hl 2v q^^ n ), we have 0 < (v,n) < 2. 
Thus, the vector v can be locally decomposed in two parts: v = ((v, n ), (v, x)) = (v n , v^) ■ 
Thus, returning to N(t) , we obtain with the error 0(t) which comes from the Jacobian 
approximation (see |J(s, 0)| in Section H]) 

N (t) = I -^T e ~^ 2dv x [ —i=e~^ 2 ( f X 2 ^Di,v n ( x ) dx ) dv n + 0(t ) 

Jr™- 1 7t 2 j o \Jtx \Jn e / 

p2 -z 2 

= / < -—=/i(dVL,2'/Dtz)dz + o(t^~). 

Jo a/7T 

If <9hl is regular, then d = n — 1 and ) = o(y/i ), which, as it was mentioned, is 

actually 0(t). The last formula that depends only on a volume of the interior Minkowski 
sausage, holds for all types of connected boundaries described in Subsection 12.21 
The formula (1771) follows from Eq. (l76p and the relation 

/u(dfl, ez) = z n ~ d fi(dQ, e) + 0(e 2(n - d) ), (79) 

which, for a fixed z and e —y + 0 , is evident for the regular case and can be proved by 
approximating the fractal volume by a converging sequence of the volumes for smooth 
boundaries. For d = n — 1 in Eq. (1771) . one has j3\ = ~ 0.2769. □ 

A comparison between the asymptotic formula (1771) and a numerical solution of the 
problem (1741) — (1751) is illustrated in Fig. [4] (for a square and a prefractal domain). 
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Figure 4: Comparison between the asymptotic formula (1771) (solid line) and a FreeFem—— 
numerical solution of the problem (1741) — fl75[) (circles) for two domains: (a) the unit 
square (with Vol(<9fi) = 4) and (b) the second generation of the Minkowski fractal, 
with Vol(dfi) = 2 2 • 4. We set D + = D. = D = 1. 

6.2 General case 

Let us come back to the problem (fTj) —(TX|). 

According to Theorem [4] (Eq. ([52]) ). the heat content can be found up to the terms 

3 

either t 2 , or t (depending on values of A), by integrating over all boundary points 9 of 
the solution u hom of the homogeneous problem (1491) — (15 1 [) with constant coefficients taken 
at a boundary point (0, 6 ). Obviously, Eq. (1571) is valid only for regular boundaries. 
Let us reformulate it to allow an explicit calculation of the heat content for all types of 
boundaries mentioned in Section [71 

For this purpose, given e = 0(y/t ), e > y/AD + t , we divide dQ (which is still supposed 
to be regular) into J disjoint parts Bj (j = 1 ,..., J ) of the size h" - " 1 with 0 < 6 < e 
such that <90 = U j = ±Bj . 

For t —y +0, <5 —y 0 and thus, due to regularity of 90 on each Bjx] 0,e[ the 
local change of variables from Section |4] is a C 1 -cliffeomorphism. In addition, since a 
continuously depends on A (see Theorem [7]), u hom , considered as a function of 6 , by the 
continuity of A, is continuous on 6. Therefore, by the mean value theorem and due to 
the positivity of | J(s, 9) \, we deduce that for all j = 1,... J there exists 9q G Bj such 
that 


d6 / ds(l-u hom (s,M))|JM)| = / ds{l-u hom (s,9lt)) / d6\J(s,9)\. 


J Bj J [0,e] 

From Eq. ([69]) , Eq. (1571) becomes 

j 


'[ 0 , 6 ] 


IB, 


' Bj 


N (t) - V / ds(l-u hom {s : 9lt)) [ d6\J(s ,1 
0,4 J 1 

0(t fi(dQ, Vt)), 0 < A < 00 

0{y/t yu(9h2, Vt)), A = cxo 


Hence we prove 
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Theorem 5 The heat content for the solution of the problem can be explicitly 

found for all types of boundaries dfl (a connected boundary of a compact domain described 
in Subsection \2.2\) using the following expressions: 

1. for A < oo on dFt: 


N{t) — 


2 y/t 


y/HfV ol(<9fl) L‘ 


/i(<9f2, \JAD + t) j da\(a) / d zf(a,z,t) 


ion 


d zfi(dQ, \JAD + t{z — 1)) / daX(a)f(a,z,t) 


'an 


— / d zfi(dn, y/ AD + tz) / daX(a)f(a,z,t) 
Jo Jan 


+ 0(tp,(d£l, y/t)), (80) 


where a = J — H- and 

VO! y/D+ 


f(a, z, t ) = exp (2X(a)ay/iz + X(a) 2 a 2 t) Erfc(z + X(a)ay/t). 


(81) 


2. for X = 00 on dfl: 


N(t) — 


2y/DZ 


„ 2 ,2 

^ e -z 


y/D- + a /D + J 0 y/n 


p(dQ, yj AD + tz)dz + 0{yft p(dQ, Vt)). (82) 


Formulas JfEW and H82)) can be approximated by 
1. for X < 00 on dfl: 

2 y/t p(dQ, y/A D + t) 


N(t) — 


^DfVol(dn) 


.Jan 


daX(a) / dzf{a,z,t) 


d z(z - 1) 


n—d 


IdQ 


dcrX(a)f(a, z , t) 


dzz 


n—d 


fan 


daX(a)f(a, z, t) 


+ o(\Jt p(dQ, Vty 


2. for X = 00 on dFt: 

Nit) — 


2 sjDF Pn-d 


VDZ+VDP 

where /3 X was defined in Eq. 178 1) . 

Proof. Using Eqs. (T56li (16CT1) . N(t) becomes 

j 


p{dQ, yj AD + t) + 0(p(dQ, y/tf 


N(t) - M.dn, ( ) + Y j 

,=i Jm 2 


(83) 


(84) 


dsids (h+(s, si, t) — /+°(s, si, t)) / d9\J(s,Q)\ 


'Bn 


0(t p(dn, y/AD + t)), 0 < A < 00 

0(y/t p(dQ, y/AD+t)), X = 00 
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Let us calculate it explicitly. We start with the part 


el 


Nhj(t) = / dsids /i_|°(s, Si, t) / dd\J(s, 


'[ 0 >£ 


' Bj 


Changing variables as in the proof of Theorem HI Nhj(t ) becomes 


Nhj(t) = d9 

J Bj J 1 

+a(A, 0, 9q) [ d6 
Jb :i . 

Therefore, we obtain 


7T 


7r 


Nhj{t) = / 


'B, 


'[ 0 ,e 


l[o, e ](s)l[ 0 , e ](s - v / 4T> + tz)|J(s,6 l )|dsJ dz 
f l[ 0i£ ](s)l[ 0 , e ](—s + ^4D + tz)\J(s,9)\ds] dz. 


| J(s, 0)|ds 


7r 


'[ 0,6 


+a(A, 0, 0^) 


(l[0, e ](3) -l [O|e]+v ^^(s))|J(s,0)|ds j dz 


l[ 0 ,e](s)l [ _ e|0]+ ^i z (s)k(s^)|ds J dz 


7r 


Applying formula (1751) with v = cTTA/f z (see also Subsection 16.11) , we find 


Nhj(t) = n(Bj, e) — (1 — a(A, 0, 9 J 0 )) 

Thus, for A < oo Nhj(t) = n(Bj,e) since a = 1. 
We treat the second part in the same way, 

el, 


\JAD + tz)dz. 


N fj( t ) = ~ dsids f+(s, si,t) / d9\J(s,9)\, 


'[o,< 


' B, 


which is equal to zero for A = oo . For f(9 q, z, t) from Eq. (j8ip . we find that 

N fj(t) = - 2A( ^ [ dsdzl[o te] (s)l[o te] (-s + 2y/D^iz)f(z,t) I \J(s,9)\d9 
V L> + Jr 2 J Bj 


2A(^)y7 
VD^ 

+ [ d zf(0 J o ,z,t) 
Jo 

2\{0i)Vi 


dzf(9 J 0 ,z,t ) 

•*yj£D-\-tZ 



a/ 4fly 

7(,-1)v^ 
I J(s, 0)|dsd0 



| J(s, 0)|dsd0 


Bj JO 




K B j, V^ D +t) / /(0g,*,t)dz 


-y ^AD + t(z - l))dz + j f(9 J 0 , z, t)fi(Bj, yj 4£>+tz)dz 

Putting two results together, we obtain the following approximations for iV(i) : 
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1. for A < oo on dfl: 


N(t) = X>( B„ s/if H»i,z,t)Az 

,-i ' JU ' 3 1 

2 A 


U 

~[ f%,z,t)n{Bj, V4D + tz)dz 


f%, Z, t)n(Bj , ^/4D + t(z - l))ck 
+ 0(t /i(<912, \7)), 


2 . for A = oo on dQ: 


m) = x; 


J «2 _*2 

rZ e z 


VD- + \/-D+ . =1 Jo 
0(y/t fi(dQ, Vi)). 


n(Bj, V^D + t z)dz 


(85) 


( 86 ) 


It means that if the formulas for n(Bj,5) are known, we get the approximation of N(t) 
up to terms of the order of f” 2 +2 for A < oo, and of the order of f 1+ 2 for A = oo . 
Moreover, this approximation, depending only on the volume of <9h2, holds for all types 
of boundaries, even fractals (see Subsection [2^_2J and p. 378 of Ref. EJ for a similar 
conclusion). 

Let us now change the sum over j with the integral over z and make J —y +oo : 


j 

lim 'S'' C(z,t,6i)fi(Bj, V^D + tz) 

J —>+oo ^' 

3 = 1 



C(z, t, a) dist(cr, VAD + tz)da, 


where da is understood in the sense of the Hausdorff measure (d -measure) defined on 
<9h2. Thus, again with the help of the mean value theorem, we have 

J C(z, t, a) dist (a, yjAD + tz)da = f C(z,t,a)da, 

from which Eqs. (15U1) and (1521) follow. We use Eq. (17U1) to obtain formulas (1551) and (155j) . 

□ 


7 Regular case 

In the case of a regular boundary we provide the asymptotic expansion of the heat content 
up to the third-order term. 

In this case, we can approximate the solution of the system ()29]) - (l33]) by the solution 
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v of the following problem (instead of (135]) - (138|) . as previously) 


9_ 

dt 

d_ 

dt 


u + - D H 


U- — D- 


cP_ 

ds 2 

ds 2 


n —1 


d 2 


72 — 1 


^ 7^2 u + + D + J2ki(0 o) 


2=1 

72—1 


de\ 

d‘ 


2=1 
72 — 1 


£ ap I “-+ D -I> ((,o) 


2=1 


2=1 


du + 

ds 

du _ 
9s 


M+|f=0 — 1, W_|t=o — 0, 
dzt— 

D— — ls=— 0 ^+)|s=0; 


D , 


' 9s |s= 
du+ 


9s 


|s=+o — D_ 


du _ 
9s 


I s=—o- 


0, 0 < s < e (87) 

0 , —e < s < 0 (88) 

(89) 


(90) 


In this approximation the remainder terms of the system (|29)) - (l33p contain only the 
coefficients of the order \[t (to compare with (170)1 ): 


72—1 

-R(si, &o) = s±D + k 2 {0 o) + 0(s 2 ), 

i =1 


that gives 


72—1 


R(si, 0 o ) = (s T 2 y/D+(t - t)z 1 )D + Y, k-(^) + 0(t) = 0(y/i). 


i =1 


The basis of the parametrix is the Green function given by (see Section m 


h+(si, s 2 , t) = 


exp 


(si - s 2 - tD + 'y(9 0 )y 


+a( A, 0, Qq) exp 


y/4nD + t V 1 V AD+t 

(si + s 2 - tD + 'y(9 0 )) 2 ' 


AD+t 


(91) 


(92) 


/+(si,s 2 ,t) = b{ A, 0,6» 0 ) 


KOo 

D+ 


•exp ( A( ^ (si + s 2 - tD+7(0 o )) + X(0 o ) 2 a 2 t 




■ Erfc 


si + s 2 - tD+^(6 0 ) 

— urn — + xie ° )a ^) ’ 


where 7(0 O ) = Yli=i ki(0 o ). The estimate (IQOll becomes 

N(t)-f do 0 [ d S (l-fi^Mo,i)WMo)l = { ° n fJ: (!'\ <oc 
Jdfl J[ 0,e] l b A — OO 


(93) 


(94) 


Consequently, for the regular case we have 
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Theorem 6 Let Q be a compact domain of M n with a connected boundary dfl € C' 00 (M n ) . 
Then for A = oo we have 


N(t ) = 2 


1 - e' 4 


0F yj D + + yj D_ 
In the case of 0 < A < oo, we have 


Vol(<9f2) Vt + 0(C). 


(95) 


N(t) = AC 0 t [ A (a)da - 


'an 


— y/D + (n — 1) / A (a)H(a)da 
Jan 

where H is the mean curvature, and 


VDf v/zr 

+ o(t 2 ), 


X 2 (a)da 


'on 


C 0 = 1 + -erf(l) - - erf(2) + 2= ( - —A I “ 0-2218, 


(96) 

(97) 

(98) 


Proof. Let us consider the case A = oo . Using the Green function given in Eqs. (1921) 
and (1931) . we obtain 


_ 7r \e 

1 5e~ 4 — 4e _1 

C\ = —= — 6 H-—-5 erf(l) + 11 erf(2) ~ 0.5207. 

< 7r Ji r 


p(Bj, \JAD + tz + tD + ^(6 J 0 ))dz + 0(C), (99) 


' e -r 

= / 

1=1 J « 

2 a / Jf)_ 

where (3 = f]j + /g p ■ I 11 ^q. the remainder term also contains the integrals 


f C 


J + — 


(e j } dz . From Eq. (1991) we End 
’ -2 


jv(*) = E 0 

i=i J 

+0(P). 

Therefore, we have 




r* 2 yjixftz +7 (ay) -D+* 


71 JO 


ds(l-s(n-!)#(£)) 


N(t) = Vt ( 2Cy/Dl^\o\(B j 


3 = 1 


T(n- 1) ( 

vi=i 


Vol (B j )H(x j )~ / tf(cx)da 




+ 0(C), (100) 


c = 

e= 


1 — e 


-4 





V^+V^’ 

>" 4 

— erf(2)^ 

D+y/DZ 

, P 
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In addition, for all o G Bj , the distance between Xj (which also belongs to Bj ) and a 
goes to 0 as J —> +oo. Thus, since 


H(xj) - H(a )| < H'(a)\xj - a\ < CVol(Bj), 


we have 


j 

lim > 

J-^+o 

3 =1 


Vol (Bj)H(xj) 


H(cr)da 


0 . 


Hence, from Eq. (HOOD we obtain Eq. (195]1 . 

The case 0 < A < oo can be treated in the similar way using in Eq. (1851) the expansion 
of the f(a,t,z): 


/(< 7, t, z ) = exp [2X(a)aVtz + A 2 (cr)a 2 (cr)t) Erfc(^ + X(a)aVt) 

1 


= Erfc(z) — 2X(a)aVt ( —j= e 2 — zErfc(z 

' 7T 


0 (f)- 


□ 
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A Definitions of Besov spaces on fractals 


Let us define the Besov space B^’ 2 (dfl) on a d-set dtt (see Ref. [21] p.135 and 
Ref. [22]). 

There are many equivalent definitions [23] 3Tj of Besov spaces. To give one of them, 
we introduce r'll l23 j a net A f with mesh 2~ u , v € N, i.e. a division of M n into half-open 
non-overlapping cubes W with edges of length 2~ v , obtained by intersecting M n with 
hyperplanes orthogonal to the axes. In addition, we denote by Vk(J\f) the set of functions 
which on each cube W in the net Af coincide with a polynomial of degree at most k. 

Definition 4 (Besov space B^ ,q (T), f3 > 0, see Ref. [21)/) Let T be a closed subset 
of M n which is a d -set preserving Markov’s inequality for 0 < d < n and let rrid be a 
fixed d -measure on T . We say that f e B p f q (Y), (3 > 0 , 1 < p,q < +oo , if f E l/imf) 
and there is a sequence B = (B u )6 l q such that for every net Af with mesh 2~ v , 
v e N there exists a function s(Af) € (Af) (by \(3] is denoted the integer part of (3) 
satisfying 

11/- o(Af)||„ (mj) < 2-^B v . 

The norm of f in B p f q (T) is given by the formula 

\\f\\B p p’ q (r) = ||/||ip(m d ) + inf \\B\\eq, 

where the infimum is over all such sequences B . 
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The dual Besov space (B 2,2 (dfl))' = B^(dfl) is introduced in Ref. |23| . To give the 
definition of the Besov space B_g{dil) we need to define the atoms: 

Definition 5 (Atom f 23f) Let /3 > 0, 1 < p < oo ; and let W with W fl T 0 be 
a cube with edge length 2~ u , is E N. A function a = aw G L v {md ) is a (—/3,p) -atom 
associated with W if 

1. suppa C 2W, where 2W is the cube obtained by expanding W twice from its 
center, 

2. f x 7 a(x)dm d = 0 for I'y | < [/3\ if is > 0, 

Let Afu ( T) = {W E M u | W fl r 7^ 0} with the notation M v of the net with mesh 2~ v 
such that the origin is a corner of some cube in the net. Then we can define the Besov 
space with a negative parameter — {3 , B^^dLl), which is actually [23] the dual Besov 
space of B 2 ’ 2 (dn) : 

Definition 6 (Besov space B r f r ((Y), (3 > 0, see Ref. [21$) The space B^(T), [3 > 
0, 1 < p, q < oo consists of functions f E which are given by 

e V(R n ) (f,(j>) = X] Sw [ aw4>dm,d, 

^gn we M v (r) 

where aw are (—/ 3,p ) -atoms and sw are numbers such that S = (S u )u<m G £ q and S v 
is defined by 

Su = | ^ | Sw 

\weA4(r) 

The norm of f is defined by 



\\f\\B™(r) — inf IlSlk 

where the infimum is taken over all possible atomic decompositions of f : 

f = ^ ^2 s waw- 

ueN WeA4(r) 

B Explicit computations for half space problem with 
constant coefficients 
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B.l Case A = oc 

The Green function of the one-dimensional problem (T^JT) — (T5TT) with A = oo and s6l 
was treated in Ref. EHU and it is given by 

r(s, Si, t) i{ S >o ! S i> 0 }r ++ (s, si, t) t i{s<o, si> 0 }r —hi 1 ®) ®i> t') 


with 


r++(s,si,t) = 


r_+(s,si,t) = b 


y/AnD+t 
1 


exp 


4 D:t 


+ A exp 


(s + ai) s 

4 D+t 


\/vr D + t 


exp 


/ 


V 


, d_ 

s ~ Si \/d: 
4 D_t 


! \ 


/ 


, ( 101 ) 
( 102 ) 


where A = and B = _ _ . 

\/ D — yjD+ + yJD- 

Let us use this result to hnd the Green function T re9 (s, Si, t) of the following one¬ 
dimensional problem 


d d 2 d 

—u+ - D + —u+ + D + k—u + = 0, s > 0 

d d 2 d 

—u_ - £>_—«_ + D_k—u_ = 0, s < 0 

at os z os 

u + \t=o = 1 , U-\t=o = 0 , 

I i n 1 I _ n A i 

s=+0 s =—0 5 s=+0 J-)— ^ U— s=— 0 i 

OS 


ds 


The constant coefficient problem 


d d 2 d 

—u - D——u + Dk—u = 0, s G M, 

at os z os 

u\ t= o = u 0l 

has the Green function of the form 

1 _ {s-y-tDkf 

h{a - s " t) = 7m e • 

that means that the change of variables s — tDk = X transforms (11071) to 

d n d 2 . 

Wt u ~ D SP“ = 0 


with the Green function 


1 (X-Y) 2 

A '»™ = 4 


(103) 

(104) 

(105) 

( 106 ) 


(107) 

(108) 
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In addition[6], we know (see (I101IMI102D I the Green function for the constant coefficient 
problem 

d d 2 

di u+ ~ D+ dx* u+ = 0 x>0, 
d d 2 

m u - ~ D -dT* u - = 0 x<0 ' 

u + \t=o = 1 , u~\t=o — 0 , 

I | n 9 I - n 9 i 

M+|x=+0 ~ U-\x=- 0 , ^ + ~QX U + |x=+0 — lJ-T^7U-\x=-0- 

Consequently, we perform the following change of variables in Eqs. (ll03IMI106lk 
X = l s>0 (s) (s + tD + k) + l s<0 (s) \j=j~ s _ tD - k 


and obtain for z = j^X that 


Thus, 


d d 2 

dt u+ ~ D+ dx* u+ = 0 x>w + k > 

d d 2 

z<tD + k, 

u + \t=o = 1, U— |t=o = 0, 

d 


d 


U + \x=+tD + k ~ U-\ z =-tD+k, D + -^U + \x=+tD + k ~ D-—U-\ z= - t D + k- 


dz 


r+^(s,si,t) 


+Aexp 


exp 


(s — Si — tD+k)'* 


4 D+t 

(s + Si — tD + k ) 2 


4TU 


for s, Si > 0 , and 


T 7 i e f(s, si,t) = 


y/irD+t 


B exp 


{[D + /DJ\s - s\yjD + /+ tD + kf 
4 D_t 


for s < 0, Si > 0. Now, to obtain the Green function of the multidimensional problem 

n—1 


a_ 

at 

d_ 

dt 


u+ - D h 


— _D_ 


a^ 

as 2 

a^_ 

as 2 


a^ 

3«, 2 

k!_ 

so. 2 


E 

i— 1 
n—1 

E 

i— 1 

M+|t=0 = lj M_|i = o = 0, 

a a 

^+|s=+0 M-| s =-0, -D+|s=+0 -^~U— |s=_ 0; 


a 

M + + D + k—u + = 0, s > 0, Qi G M, 
Os 

a 

u_ + D-k—U- = 0, s < 0, Qi G M 
os 
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we apply the Fourier transform in s, variables and, due to the boundary conditions 
depending only on s , we obtain that G±+ , the Fourier transform of the Green function 
G± + , can be found by the formula 

G±+(s,Si,f,t) = e _D±l 5 | 2 t r ±+ (s,Si,f), 

where r± + (s, si, t) is the Green function of the corresponding one-dimensional problem. 
This implies 


+Alexp 


(4:7TD + t) 2 
(s + Si — tD + k) 2 


4 D+t 


exp 

exp 


(s — si — tD+kf 
4 D+t 
\9-9^' 


AD+t 


for s, si > 0, 9 , 9\ G 


nn—1 


rh+(s,0,si,0i,t) = 


B 


(AirD_t) n 2 y/nD + t 


exp 


|r s _ s i + tD + k) 2 
4 D-t 


• exp 


|0-0i| 

ADB 


for s < 0, si > 0, 9 , 9i G 


nn—1 


For T and Y rea we also have Varadhan’s bou n ds[30 ] for s ^ si 


lim t lnTr , (s, Si, t) 
o+ 


lim fin 

t —^0-1- 


r+f(s,si,f) 


lim f In r_ + (s, si, f) = lim t In ry+fs, Sj, f) 

t->0+ t-y 0+ 


d(s,s i) 2 
AD + 



2 


4L>_ 


where d(s,s i) is the Riemannian distance between s and Si, which is equal here to the 
Euclidean distance, since D + and D are constant in Q_ and Q_ respectively. 


B.2 Case 0 < A < oc 

Let us consider the one-dimensional problem f|T9j) -([5T ]) with A = A(0o) and s G R. The 
associated problem for the heat kernel is then given by 

(d t - D±d 2 s ) G(s,si,f) = 0, 

G\ t =o = <5(s, si) for s > 0, 

G(—0, si, f) = A(G(—0, si, f) — G(+0, si, £)), (109) 

D_j_—G(+0, Si, f) = G(—0, si, f). (HO) 

We search the explicit solution of the problem[h] with 

G(s, Si,f) = | Q+ + 
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s < 0, Si > 0 
s > 0, si > 0 












We seek for G_+ and G++ in terms of free heat kernel K(s, s 1; D±t) (see Eq. (l55j) l and 
single layer heat potentials for Si > 0 : 


G ++ (s,si,t) = K(s,s u D + t) + D + K(s,Q,D + (t - r))a+(si, r)dr (s > 0), 

Jo 

GL+(s, si,t) — D_ / K(s, 0, D-(t — r))a_(si, r)dr (s < 0), 


where «±(si,t) are unknown densities to be determined. Considering the boundary 
conditions (1 109jl — ill 10|) and the jumps of the hrst derivatives of G± + , 

d Id 

^^++|s=+o — — 2 a+ ^ Sl ’C + ^-^(0) s i) D+t)i 

d 1 

—G__|_| s= _o = --a-{si,t), 


we obtain two relations 


d 


D_a_(si,t) = -D + a + (si,t) + 2D + —K(0,s 1 ,D + t), 

ds 


D_ot-(si,t) — 2XK(0, si, D + t) + A 


VD+ f* a+(si,r) 


Jo y/t-T 


dr 


-A 


yjDJJ f* Q-(si,r) 

7T Jo yjt-r 


dr. 


Following the method from Ref. [ 6 ], we solve the system corresponding to a_(si,£) 
and a + (si, t ): 


D_a-(si,t) + D+a+(si,t) = 2D+JJ-K (0, s u D+t), 

ds 


-ZA_Q!_(s^, t) ~\~ 


X yfDZ 


7r 


1 + W ^7 


Q_(si,r) 

y/t — T 


dr = 4AA"(0, si, D + t). 


We obtain therefore the Abel integral equation of the second kind for a_(si,i) 

a_(si,t) + 7 


f ‘ t a_(si,r) 4A 


- dr = — AT(0,si, D + t), 

y/t — T A>_ 


where 7 = 


y/^ D - 


1 + 


Consequently, 


4A 


«-(si ,t) = —K(0,S!,D + t) -1-Q- j 
+vr 7 2 -^ /V' 2(i - r) (k(0, si, D + t) — 7 


4A [* K{Q,s u D + t) 


dr 


y/t — T 

' T K^XsiGJ+s) 

y/r — s 


ds dr. 
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Using the Laplace transform yields, after simplifications: 

1 ( ( (s-Si) 2 


G++(s,s\,t ) — 


y/Air D + t 
A 


D. 


exp 


where a = 


1 + 1 


/D_ 


'D 4 


exp 


Act 

7K 


AD + t 


+ exp 


(s + si) + A a t Erfc 


(a + ai) s 
4 D+t 
f s + Si 

VV^ 




+ 


A aVt ) , 




By the same way, 


G_ + (s,Si,f) — 


A 


f Act ( 


v^nx exp \7Vz {- s + s '\ll5: 


+ Va‘t 


— S + Si \/j— 

Erfc | -r— + + Xay/t 


2 y/Dlt 


We see that the Green function G ++ for A = 0 becomes the Green function of the 
problem with the Neumann boundary conditions and in this case N(t) = 0, as w_ = 0. 
This property, N(t) = 0, can be also directly found using the Green function. 

In M n for x = (s, 9) and y = (si, Of) gRx M n_1 we have 

G++(s,9, = G++(s, 

G —h( s ) s i 5 OiG)«. n — G~ + (s, si,t)^K(9,9i, D_t)^n~i. 

Therefore in M n for Varadhan’s bounds with x ^ y we have 

d(x,y) 2 


lim + 1 In G ++ (x, y, t) R n = 


lim thiG- + (x,y,t)-^n = 


4 £>_, 
d(s, Sii 


+ d(0, 9\ 


AD_ 


Remark 3 Applying this framework to the same system but with the transmittal boundary 
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condition for 0 < A < oo, we obtain 


G++(s, 9, s\,9i,t) — 


exp 


(47r D + t) 2 
(s + Si — tD + k) 2 


exp 


(s — si — tD + k)‘ 


4 D+t 


ADM 
A 


exp 


9 ( 6 , 6 ,)' 

ADA 


(AnD+t)™? D + 


exp 


/ Act 

WK 


(s + Si — tD + k ) + A 2 a 2 t 


• Erfc [ S + Sl ^!F +k - + XaVt) exp 


2y/D+t 
G — |_(s, 6 , si, 9,,t) = 


9(6,6,)'- 

ADA 


A 


(47riA_t)"2 y/D_D + 


■ exp | —f== ( + Sl \lfj- + tD+k ) + A 2 a 2 t 




■ Erfc 


— TT^S + Si \!+ tD+k 


2y/DZt 


A aVt 


exp 


9(6,6,)' 

ADD 


We also notice that for a fixed t > 0 for A —> +oo we obtain 

G ++ (s, s,, t) —> T ++ (s, si, t) and G_ + (s, s,,t) —)■ T_ + (s, s,,t). 
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